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Density Matrix Renormalization Group (DMRG) calculations on 4-leg t — J and Hubbard ladders 
have found a phase exhibiting "stripes" at intermediate doping. Such behavior can be viewed 
as generalized Friedel oscillations, with wavelength equal to the inverse hole density, induced by 
the open boundary conditions. So far, this phase has not been understood using the conventional 
weak coupling bosonization approach. Based on studies from a general bosonization proof, finite size 
■ spectrum, an improved analysis of weak coupling renormalization group equations and the decoupled 

' 2-leg ladders limit, we here find new types of phases of 4-leg ladders which exhibit "stripes". They 

also inevitably exhibit "bipairing", meaning that there is a gap to add 1 or 2 electrons (but not 
C , 4) and that both single electron and electron pair correlation functions decay exponentially while 

, ^ ' correlation functions of charge 4 operators exhibit power-law decay. Whether or not bipairing occurs 

' in the stripe phase found in DMRG is an important open question. 

in : 

I. INTRODUCTION 

"^i . . . 

O , There is clear experimental evidence for "static stripes" i.e. charge density waves (CDW's), in some cuprate 
^ ' superconductors at commensurate doping (where the superconductivity is suppressed) [lIQ. It has been suggested 

] that "fluctuating stripes" may occur at incommensurate doping and may be responsible for superconductivity. It has 
S^' also been suggested that stripes only occur in models where the long range Coulomb interactions are kept 3 • the 
^ ■ other hand, numerical evidence for stripes has been found in Hubbard and t-J ladders, which contain only short range 
, t , interactions An understanding of the occurrence of stripes in these models, whether or not long-range 

' Coulomb interactions are required for their existence and their connection with superconductivity are important open 
^ , questions. 

There are many experimental realizations of ladder systems (lllfl^ islfl^ . Recently, stripe has been observed by 



a resonant X-ray scattering technique in a 2-leg ladder compound Il5l|. Numerical evidence for stripes comes from 



Q • density matrix renormalization group (DMRG) work [HsllGllTllsllglllcil. Since this method, as originally formulated, 
^ ^ is intrinsically one-dimensional (ID), the results have been presented for "ladders" i.e. finite systems in which the 
I— number of rungs is considerably greater than the number of legs. (For instance, results exhibiting stripes have 
been presented for systems of size 6 x 21 [lo|.) In the limit where the length of the ladders (number of rungs) is 
^ , much larger than their width (number of legs) they become ID systems and a corresponding arsenal of field theory 
J — ' methods, such as bosonization, can be applied. Combining DMRG results with field theory methods to extrapolate to 
', the limit of infinitely long ladders is an important step towards reaching the two dimensional (2D) limit. Of course, 
t — ■ ■ an extrapolation in the number of legs must finally be taken. 

■ DMRG works much more efficiently with open boundary conditions (OBC) and most work on ladders has used 
\^ OBC in the leg direction. Such boundary conditions can induce "generalized Friedel oscillations" , meaning oscillations 
in the electron density which decay away from the boundary with a non-trivial power law and oscillate with an 
, incommensurate wave- vector, often related to the hole density p^ . While sometimes regarded as an unphysical 
• nuisance, we regard OBC as a useful diagnostic tool. According to bosonization results, the density-density correlation 
J> \ function for an infinite length ladder decays with twice the exponent, the same wave- vector and the square of the 
" ■ amplitude, governing the Friedel oscillations. Thus the Friedel oscillations are giving information about correlations 
rxj ^ in the infinite system. Furthermore, in the 2D limit of an infinite number of legs, these Friedel oscillations could turn 
^ I into a static incommensurate CDW or else fluctuating stripes. 

■ ■ ' While a true long-range CDW is possible at commensurate filling even in ID, bosonization/field theory methods 
suggest that it is not possible at incommensurate filling in ID, giving way instead to boundary induced Friedel 
oscillations [TgHTtI . (A long-range CDW at p/q filling, where p and q are integers, becomes increasingly suppressed as q 
increases.) This assertion is related to Coleman's or Mermin- Wager's theorem about the impossibility of spontaneous 
breaking of continuous symmetries (and the impossibility of the existence of the corresponding Goldstone modes) 
in Lorentz invariant ID systems even at zero temperature (or correspondingly 2D classical systems at any finite 
temperature). A true long-range incommensurate CDW leads to spontaneous breaking of translational symmetry 
by any integer number of lattice spacings, no matter how large, whereas a commensurate CDW only breaks a finite 
dimensional symmetry. Taking the low energy continuum limit, translational symmetry at incommensurate filling is 
promoted to a true continuous U{1) symmetry and Coleman's theorem apparently applies. 

"Stripes" or boundary-induced Friedel oscillations, have been observed in the t-J model, related to the U —f oo 
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limit of the Hubbard model, on 2-leg ladders [l^. These Friedel oscillations, at a; ^ 1, are of the form: 
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a=l 




(1.1) 



Here the average electron density is. 



n = N,/{2L), 



(1.2) 



where Ng is the total number of electrons, i^+p is the Luttinger parameter for the charge boson (p) which is the sum 
(+) of the 2 charge bosons corresponding to the two bands in a weak coupling analysis. A and a are constants. Note 
that we may replace the oscillation wave- vector by: 



where 5 = 1 — n, is the hole density, measured from half-filling, since x is always integer. "Snapshots" of typical 
configurations within the DMRG calculations show pairs of nearby holes, one from each leg, well-separated from other 
pairs. An appealing picture is that the holes are pairing into bosons, 1 hole from each leg. (Here we refer to bosons 
with a conserved particle number, such as atoms, not the bosons arising from bosonization.) This phase is of CISO 
type, indicating that only 1 gapless charge boson survives and zero gapless spin bosons, out of the 2 charge and 2 
spin bosons introduced in bosonizing the 2 leg ladder. This phase exhibits exponential decay for the single electron 
Green's function but power law decay for the electron pair Green's function. One may approximately map the 2-leg 
fermionic ladder into a (single leg) bosonic chain. The standard superfluid phase of this boson model exhibits Friedel 
oscillations at wave- vector 2TrS where 5 = N^/L, is the number of bosons (i.e. hole pairs) per unit length and the 
number of bosons is: Nb = Nfi/2, were Nh is the number of holes. These Friedel oscillations in the bosonic model 
just correspond to a sort of quasi-solid behavior. We can think of the bosons as almost forming a solid near the 
boundary with a uniform spacing between all nearest neighbor bosons. This, of course, coexists with quasi-superfluid 
behavior since the phase correlations also decay with a power law. Thus it could be called quasi-supcrsolid behavior 
and is typical of many ID systems. In general the density profile will contain other Fourier modes besides the 27m 
mode kept in Eq. (|1.1|) . This 2TTn mode dominates in the sense that it has both the smallest wave- vector and also 
the smallest power law decay exponent. Note that this behavior is quite different than what we might expect in a 
C2S0 phase, for example, or which occurs at zero interaction strength in the C2S2 phase. Then we expect density 
oscillations at wave-vector 2kp^e and also 2fci?^o where k p^./^ are the Fermi wave- vectors for the even and odd bands. 
(They are even or odd under the parity transformation that interchanges the two legs.) In the CISO phase that is 
observed in 2-leg ladders, the oscillation wave-vector can be written as 21m = 2(kp^e + kp^o) = '^kp where kp is the 
average Fermi wave-vector. 

DMRG works on the doped 4-leg t-J model exhibited two phases, both of which appear to have a spin gap [rllsl fislfigl . 
At low doping, the dominant Friedel oscillation wave-vector appears to be 47rn, where n, the average electron density 
is now: 



Above a critical doping Sc-, corresponding to 5c ~ 1/8, for J = 0.35t and 0.5i , the oscillation wave-vector changes to 
27rn. DMRG "snapshots" of typical configurations suggest well separated pair holes in the lower density phase but 
4-hole clusters (1 hole on each leg) in the higher density phase. This is consistent with the 27rn Friedel oscillation 
wave- vector since the average separation along the ladder of the equally spaced 4-hole clusters would be 1/n. This 
higher density phase with 2?™ oscillation wave-vector has been referred to as a stripe phase. 

In the standard weak coupling approach we assume that all the interactions are small compared to the hopping. 
Thus we first solve for the band structure of the non-interacting model and then take the continuum limit of the 
interacting model, yielding right and left moving fermions from each band. These continuum limit fcrmions are then 
bosonizcd. Letting kpi be the Fermi wave-vector of the 4 bands {i = 1,2,3 or 4) the number of electrons in each band 
(summing over both spins) is: 



27rn 



-27r(5, 



(1.3) 



nEE iVe/(4L). 



(1.4) 



Nl = L{2kp,l^), (1=1,2,3,4). 



(1.5) 



Thus we see that the electron density is: 
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(1.6) 



i=l 
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27m = ^fcF» = 4feF. (1.7) 

1=1 

Thus the stripe phase again corresponds to Friedel osciUations at a wave- vector of Akp. Using the equivalence of 2ttti 
with —2tt6, the physical picture of the stripe phase is well-separated clusters of 4 holes (one on each leg). If such 
4-hole clusters are equally spaced we obtain a Friedel oscillation wave- vector of 27rn, since the average separation along 
the legs of the clusters is 1/n. On the other hand, the lower density phase shows no evidence for Friedel oscillations 
at wave-vector 2Tm and instead Aim {8kp) oscillations appear to dominate. (We note that, when interactions arc 
included, we expect the Fermi wave-vectors to be rcnormalized or lose their significance entirely. However, the sum of 
all Fermi wave-vectors, Akp = 2Tm, is known to still be a meaningful and "unrenormalized" wave-vector even in the 
presence of interactions. This follows from the ID version of Luttingcr's theorem, proven in Ref. [20|). This stripe 
phase appears to have a spin gap. Limited DMRG results have been presented on the decay of the pair correlation 
function. Pairing correlations appear to go through a maximum, as a function of doping, at a somewhat higher doping 
than 6c where their behavior appears consistent with power-law decay 001 . 

A simpler picture of the stripe phase of 4 leg ladders was obtained by mapping each hole pair on the upper 2 legs 
into a boson moving on a single chain, and likewise mapping each hole pair on the lower 2 legs onto a boson moving 
on a different single chain. Thus the 4-leg fermionic ladder is mapped into a 2-leg bosonic ladder p^ . 

Bosonization analysis of the 4-leg ladder it is not an easy task. A total of 8 left and right-moving bosonic fields 
must be introduced, one for each spin component in each band (or rung). In a general phase any of these boson fields 
or their duals might be pinned leading to 3^ possible phases! Taking into account that the points where the fields are 
pinned are of physical significance, leads to even higher estimates of the number of possible phases. Additional phases 
occur if it is assumed that some of the bands are completely filled (or empty). Keeping only Lorentz invariant, non- 
Umklapp interactions in the continuum limit there arc 32 different interactions terms and hundreds of non-zero terms 
in the RG /3-functions at quadratic order. Analysis of these RG equations generally indicates that some couplings 
flow to the strong coupling region where the equations break down. The bosonization analysis of [2]| predicted more 
than 10 different phases for the 4-lcg ladder with open boundary conditions in the leg direction, as hoppings and 
hole doping are varied. However, only two of these phases have a spin gap and neither of them has the right Friedel 
oscillation wave- vector to describe the stripe phase seen in DMRG work. 

An alternative bosonization approach was introduced in (2^ . This was based on the limit of low doping with small 
but finite on-site interaction U . Although details were sketchy, it seemed to correctly reproduce the two phases found 
by DMRG for the 4-lcg ladder. 

The purpose of this note is to reexamine bosonization approaches to 4-leg ladders and the nature of the stripe phase 
found in DMRG. We first give a general discussion of the conventional weak-coupling bosonization approach and of 
Friedel oscillations in 4-lcg ladders. We restrict our attention to CISO phases where there are various candidates for 
the stripe phase. We sketch a proof that all possible phases with stripes inevitably have exponentially decaying pair 
correlations within this approach. However, they can have bipairing, corresponding to quartets of fermions. The finite 
size spectrum (FSS) of the stripe phase, which indicates the total charge excitations relative to the ground state can 
only be multiples of four, is also consistent with bipairing. On the other hand, a pairing phase should correspond to 
Skp density oscillations instead of stripes. Then we show that the stripe phase could arise from an improved analysis 
of the so-called "fixed ray solutions" of the renormalization group (RG) equations. This new analysis [1^ is based 
upon the potential structure of the RG equations 2^l25t and gives the CISO phases with correlations consistent with 



our general discussions based on bosonization and FSS. 

We then discuss a different bosonization approach. We start with two identical decoupled 2-leg ladders and then 
turn on the interlcg hopping t2^± and interleg interaction V2,±, which are much smaller than the other energy scales 
in the 2-leg Hamiltonian. Now the system becomes a pair of weakly coupled identical 2-leg ladders. Each 2-leg ladder 
is in the CISO phase [1^ and we also assume that t2,± and V2,± are small compared to the gaps to the other 3 modes 
in each 2-lcg ladder. The model then becomes equivalent to the continuum limit of a 2-leg ladder of bosons, rather 
than fermions, of the type studied in [23|. This model has two phases: 47r/Oo density oscillation (with superfluidity) 
and 27rpo density oscillation (with boson-pair superfluidity) that correspond to the 2 phases seen in DMRG, where 
po is the average boson density. The Friedel oscillations correspond in the two models and the stripe phase has 
boson-pairs, i.e. bipairing. Interestingly, these two phases also have the same features as those found in the improved 
weak coupling RG, where 2ttpq, 47rpoj superfluid and boson-pair superfluid in two-leg bosonic ladders correspond to 
4:kp, Skp, pairing and bipairing, respectively, in 4-leg fermoinic ladders. 

This approach is very close to that of [l3| which maps the 4-leg fermionic model onto a 2-leg bosonic model in 
a more phcnomcnological way, based on exact diagonalization and DMRG. Our results arc consistent with those of 
that paper which did not, however, point out the occurance of bipairing in the striped phase. This approach, and 
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our conclusions are also closely related to those of [24|, as we will discuss. We point out that the numerical results 
[ist on the bosonic model which is argued to correspond to the 4-leg fermionic model in its stripe phase show some 
evidence for boson pairing. However, published results on the 4-leg fermionic ladder so far do not show a gap for 
pairs as should occur in the bipairing phase. Since the ladder length in Ref. [isj was only L = 24, the problem may 
be that the pair-gap is too small compared to the finite size gap at this ladder length. More DMRG results on longer 
ladders may be required to confirm that the stripe phase indeed has bipairing. It is similarly unclear whether pair 
correlations decay exponentially, as would occur in a bipairing phase, or with a power law as in a pairing phase. More 
DMRG work should shed light on this question. 

If this bipairing assumption is correct it finally yields a remarkably simple picture of the stripe phase in 4-leg 
ladders. It is a phase in which electrons do not form pairs, but rather bipairs. While the usual pairing does not occur 
in this phase, it tends towards a more exotic form of supercondutivity based on condensed charge-4 objects. There 
are proposals in the literature that some systems like a frustrated Josephson junction cha in I28ll29l| . strongly coupled 
fermions [sot . cuprates with strong phase fluctuations 31 1 and a spin 3/2 fermionic chain |32| can have unusual pairs 
composed of 2 Cooper pairs or four-particle condensations. Experimental evidences that Cooper pairs tunnel in pairs 
across a (100)/(110) interface of two c?— wave superconductors are also reported (33ll33 |. They are very similar to but 
not exactly the same as the bipairing here. Starting with a four band system, the phase we propose is more like a 
generalized Luther- Emery phase with four- holes. 

We note that the stripes observed in 6-leg ladders (olfiot appear to have a more dynamical character, containing 
less than 6 holes per stripe. Thus they may exhibit behavior closer to that proposed in the 2D case. 

In the next section we review some features of bosonization of 4-leg ladders including Friedel oscillations and sketch 
a proof of one of our main results: within the usual bosonization framework any phase with stripes does not have 
pairing (but may have bipairing) . This result applies very generally and docs not depend on the detailed form of the 
Hamiltonian. In Sec. Ill wc study the finite size spectrum corresponding to a stripe phase. In Sec. IV we study weak 
coupling RG equations, showing that a stripe phase can arise from studying fixed ray solutions. In Sec. V we discuss 
the case of small t_L,2, Kl.2 and the analogy with a 2 leg bosonic ladder. Sec. VI contains conclusions. In Appendix 
A, we give some more details about the proof that stripes and pairing can not coexist. In Appendix B we derive in 
detail, in the weakly coupled Hubbard model, how the density operator obtains terms leading to Akp oscillations. The 
initial conditions of RG equations in terms of the bare interactions and Fermi velocities are given in Appendix C. 



II. CONTINUUM LIMIT AND BOSONIZATION 



A. Bosonization 



In this section and the next we focus on a weak coupling treatment of the Hamiltonian on a ladder with N legs and 
L rungs with open boundary conditions in both rung and leg directions. The Hamilonian is H = Hq + Hint-, with 



Q = ± 



X-l N 



L N-1 



E E *cl,a(a;)ca,a(a- + 1) + E E *-Lc1^„ (x)ca+i,a (a;) 



a; — 1 a—1 



x—1 a—1 



h.c. 



(2.1) 



and 



N L 



Hint = E E '^a.T (^) 



(2.2) 



Here Ua^a = c\ ^Ca.a- We first diagonalize the rung hopping terms in the Hamiltonian by transforming to the band 
basis, ipi,a, for i = 1, 2, 3, 4: 



where 5 is a unitary matrix: 



(2.3) 



a=l 



(2.4) 



The dispersion relation for the j^^ band is £j{k) = —2tcoak — 2t±coa(kyj), where kyj = jir/b. We may take the 
continuum limit, to study the low energy physics, by introducing right and left moving fermions fields (sometimes we 
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call them chiral fcrmions) which contain wave- vectors of the original band fermion fields near the Fermi points, 

i^U^) - i^Rio.{x) e^'^--^^ + ^L,o.{x) (2.5) 
where kpj is the Fermi wave vector for band j . Then we can bosonize these right and left fcrmions by the dictionary 

i>R/L^c. ~ 77,e'^^'^«/-° (2.6) 
where r]i are Klein factors with {rji^iy} = 25ij. The commutation relations of the boson fields are 

1 

[<PRia(,x),ipL](}{y)] = -^SijSafl, (2.7) 

1 

[ipRia{x),ipRjp{y)] = -sgn{x - y)SijSa[3, (2.8) 
— z 

[(pLiaix),(pLjpiy)] = —sgn{x - y)5ij5ai3- (2.9) 

Now we have two chiral bosons and it's more convenient to describe physics by the conventional bosonic field (j) and 
its dual field 6*, 

0iQ = 'fRia + '■PLia,&ia = V Ria " </5LiQ, (2-10) 

since the partial derivatives of 9 and are density and current, respectively. We can also separate the charge p and 
spin (T degrees of freedom by introducing the following fields, 

(t'ip = —7^{4>i-\ + (t>il)i (2.11) 



V2^ 



fit - Vii 



) (2.12) 



and similarly for 9 fields. Then Eq. (|2.ip becomes 

Ho = Y.T f ^^[(^-'^-)' + (2.13) 

where i/ = p or cr and Vi is the Fermi velocity of band i. We define new fields from the linear combinations of the 
bosons from different bands: 

</'S"^ = ^(4p,.±0,p,.), (2.14) 

and the same for 9^f^ and 9^^^ . Eq. ((2l4l) is convenient when we consider the bosonized interactions but is not 
necessary the best basis to describe the system. The basis to describe the system should be the one in which boson 
fields get pinned, in principle determined by the interactions. We know that there should be four mutually orthogonal 
charge and spin bosons for 4- leg ladders. That is to say, starting with bosons in the band basis, Eq. (|2.1ip and (|2.12p . 
guided by the interactions, we can find the proper new basis in which some linear combinations of band bosons are 
pinned. The new and band basis only differ by an orthogonal transformation. In general, the transformations of 
charge and spin fields don't have to be same. 

Due to the different Fermi velocities in Eq. (|2.13p . after changing to a new basis, there will be mixing terms 
between the derivative of boson fields. These mixing derivative terms are only up to quadratic order in the boson 
fields. Spinless fcrmions on 2-leg ladders were studied in Ref. [s^ and they found that the mixing terms only modify 
the exponents of correlation functions within the conventional bosonization analysis. We are mainly concerned about 
whether an operator is power law or exponentially decaying but not its exponent. In the later sections, we only discuss 
the mixing terms if it's necessary. 

B. Stripes and Bipairing 



In this section and Appendix A, we wish to sketch a proof of a general result that would apply to any of these phases 
that are candidates for a stripe phase, regardless of what basis for the boson fields we use. Stripes are incommensurate 
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density oscillations with the lowest wave- vector 2im (e.g. no Trn density oscillations) [7||8|]. The wave- vector 27rn is 
equivalent to Akp via ID Luttinger theorem [23|. Stripes thus can be described as the generalized Friedel oscillations 
induced by the boundary A derivation of the Akp components of density operators is given in Appendix B. At 
incommensurate filling, the ID Luttinger theorem suggests that at least one gapless charge mode always exists. Stripe 
phases appear to have a spin gap in DMRG works and we expect that phases with more than one gapless charge 
mode are generically unstable. So we restrict out attention to CISO phases. 

We then argue that any phase exhibiting stripes (i.e. 27m oscillations) cannot exhibit pairing (i.e. power law pair 
correlations and gapless pair excitations) but can exhibit bipairing (charge four operators). 

In a CISO phase, 7 out of 8 boson fields are pinned due to interactions. We will consider all possible four fermions 
interactions written in terms of right and left moving fields 4'r/l i including those that are ignored in the conventional 
weak coupling analysis. Similar to what we discussed in the previous section, after bosonizing, if the interaction is 
relevant, the boson fields in the interaction will tend to be pinned to constants. Before we proceed, we should explain 
explicitly what we mean by a boson field 9 (or (f) being pinned. Consider a vertex operator of a boson field, e**^, 
where a is a constant. If (e"'^) = const ^ 0, then we say 9 is "pinned" to some constant modulo 2-k /q. On the other 
hand, we say 9 is "unpinned" if ^e*^^) = 0. There could be three situations for 9 to be unpinned. First, 9 is gapless. 

Second, a part of it is gapless and the rest is pinned. For example, if {9ip + ^?2p)/\/2 is gapless and (03p -f 6'4p)/-\/2 is 
pinned, then ^e"'®^'') is still zero since we can replace the pinned fields by their pinned values. Thirdly, its dual field 
(j) or part of its dual field is pinned and 9 will fluctuate violently. 

There are many such phases characterized by which fields are pinned. In general, it is not appropriate to simply 
label the pinned fields as or 9i^ in terms of the band basis. The interaction terms in the Hamiltonian involve 
various linear combinations of these fields and it is generally necessary to characterize phases by first going to a 
different basis of boson fields: 

9\^ = Rp9p, = Rp$p, 9'^ = rX, = rJ„, (2.15) 

where 9p ~ {9ip,92p,93p,94p) etc and Rp and Ra are orthogonal matrices. We will refer to the basis {9[p,<j)[p) and 
{9[aT^'ia) "pinning basis". One of the 4 charge bosons in the new basis, {9[p,(j)[p) remains gapless in a CISO 

phase. All of the other 7 bosons are gapped. However, for each of these 7 bosons we must specify whether it is 9'^^^ or 
<f>[jj which is pinned. We refer to this choice of which bosons arc pinned as a "pinning pattern" . Each pinning pattern 
corresponds to a distinct phase. (Actually the number of distinct phases is even larger than this since the points, Ci^ 
at which a boson is pinned, eg. (0;^) — Ci^,, can also characterize the phase.) 

Following conventional bosonization analysis, we will assume that all CISO phases are characterized by such a 
pinning pattern. Our basic strategy is to check through all possible pinning patterns for a CISO phase and see if any 
pinning pattern can make the correlation functions of pair and Akp density operators decay with a power-law at the 
same time. Given a pinning pattern, it's easy to know whether an operator has power-law decaying correlations or 
not. Wc can replace the pinned bosons by constants inside the correlation function. The correlation function will 
decay exponentially if the vertex operator contains the dual of a pinned boson. It will decay with a power-law if the 
vertex operator only contains the gapless mode and pinned bosons. However, there are three main complications here. 
The first one is: what is the basis for the boson fields (including the gapless mode)? Even though we know the basis, 
there are still too many pinning patterns (more than 2*"). The last problem is that the wave-vectors of Akp density 
operators can be changed and actually correspond to that of stripes if some Fermi momentum renormalization occurs. 
For example, the wave-vector kpi + kF2 + 2A;f3 = 4fcF if the condition kps = kpA is satisfied. In this subsection, 
we try to solve these problems. As a matter of fact, Rp and R^ are not arbitrarily orthogonal matrices and they 
have to satisfy some constraints due to the symmetries. Then we can classify the pinning patterns according to the 
possible gapless charge modes. For each case with the different gapless charge mode, we can systematically examine 
whether pairing and stripes can coexist or not by the inner product argument which will be introduced later. When 
the pinning patterns with all the possible gapless charge modes are checked, we complete the proof. In order to know 
what field can be a gapless charge mode, we should study symmetries first. 

We now want to discuss the consequences of symmetries. We would like to introduce a specific combination of band 
bosons Eq. (|2.1ip and (|2.12p . which is relevant to symmetry. Define the total charge (or spin) field as 

Qiu = \{eiu + 92,. + 9:,, + 9i,), (2.16) 

where v = p or a and similarly for if we replace 9 by (j). 

We can somewhat restrict the possible pinning patterns by symmetry considerations. Charge conservation symme- 
try, ij) e^'^ijj all fermion fields, ijj, corresponds to the translation: 



(2.17) 
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Thus, the allowed pinned must be the linear combinations of (f)^^ , in other words, any pinned (^-^ fields must be 
orthogonal to $ip. Similarly, the subgroup of SU(2) symmetry, ipa — > e^'^'^ipa where a = + or — for spin up or down, 
respectively, corresponds to: 

4>ja 4>ja + VV^e- (2.18) 

Therefore, any pinned can only be the linear combinations of (f>'J~ in order not to violate the SU (2) symmetry. 
We also expect translation symmetry to be unbroken at incommensurate filling, as discussed in Sec. I. We see from 
Eq. (j2.5p . that translation by one site: x ^ x + I, corresponds to the symmetry: 

V'i?/Lj.->e±'^-^Vfl/Lj., (j = 1,2,3,4) (2.19) 

corresponding to: 

0,p^0,p + ^2/^kFj, (j = 1,2,3,4). (2.20) 

Consider a vertex operator of a boson field, e*"?^, where g is a constant. Since in general the factors ^/Ij-nkFj are not 
zero modulo 27r/g, this may forbid pinning of any of the Q'^^ fields. However, it may happen, due perhaps to some 
renormalization phenomenon, that 2 or more of the kpj are equal. In that case it might be possible for some of the 
B'^p bosons to be pinned. For example if kpi = kpji then O^^J could be pinned. We will allow for that possibility since 
there is no reason to exclude it in the strong coupling region. However, can never be pinned since it transforms 
under translation by one site as: 

eip ^ Gip + 4fcF/%/2^, (2.21) 

and A.kp/i'K = n is always non-zero modulo iir/q at incommensurate filling. The fact that Oip can never be pinned 
is not exactly equivalent to the statement that the gapless boson in a general CI SO phase at incommensurate filling 
has to be (Oip, ^ip). There could be exceptions if some pinned 9'^^ fields are not orthogonal to 8ip. If that happens, 
then 9ip can't even be chosen as a basis field. As we know, all the pinned 0-p must be orthogonal to <i>ip due to the 
charge conservation. However, translational symmetry doesn't demand all the pinned O'^p must be orthogonal to 8ip 
and in general whether a 9'^p is pinned or not depends on the Fermi momentum. Thus to know whether a 6[p can be 
pinned without violating translational symmetry, we have to discuss the possible Fermi momentum renormalization. 

However, even though some 9[p field is allowed to be pinned by symmetry, in practice, whether it really gets pinned 
or not depends on the interactions in the Hamiltonian. In the weak coupling treatment, the renormalizations of kpj 
are assumed not to happen and kpj are in general all different. The interactions involving 9'^p fields won't be present 
in the Hamiltonian due to the fast oscillating factors in front of them. However, in the strong coupling region, if some 
renormalization of Fermi momentum occurs such that the oscillating factor becomes a constant, then new interaction 
will appear in the Hamiltonian. For example, if 2{kpi -\- kpj) = 27r, the interaction involving a 9[p field, such as 
g-t2(feFi+fei7j)2;^t^^^^t^^_^^_^^^.^ ^ g-«\/47r(efj ±4>ij ) ^ pg^j^ prescut in the Hamiltonian since the oscillating factor 
^~i2(kFi+kFj)x bgcomes a constant at each lattice site x. In this case the shift of 9^^ under translation is 27r/\/47r, 
which becomes 2tt for \/Att91^ . Thus if this interaction is relevant, then then the field will be pinned. 

This is an important point for our proof. We will discuss whether a 9[p can be pinned based on the Fermi momentum 
renormalization that results in the new interaction involving 9^p in the Hamiltonian. The reason is that the number of 
different oscillating factors is finite and we can certainly check all of them. In the proof, we first want to know what 
could be the gapless charge mode. 0ip and $ip are always unpinned but this doesn't mean that the gapless mode 
has to be (0ip, 'I'lp). We have to know when 9[p fields not orthogonal to Qip can be pinned. We can find all such 
situations by checking through the interactions containing 9[p fields not orthogonal to Qip and then conclude what 

are the gapless charge modes. For example if 91^ , not orthogonal to Oip, is pinned, then the gapless mode won't be 
0ip though Oip is still unpinned. This can happen if 2{kpi + kpj) ~ 2tt, as we discussed above. So we just have to 
check those Fermi momentum renormalizations which will make the interactions, containing 9[p fields not orthogonal 
to Oip, appear in the Hamiltonian. 

Besides these issues related to symmetry and interactions, the renormalizations of Fermi momenta have another 
effect on our discussions. Our goal is to answer whether stripes and pairing can coexist or not. Stripes correspond to 
the 4fcF density operators. However, a general Akp density operator may correspond to the 'ikp one if some suitable 
Fermi momentum renormalization happens. For example, the wave- vector 2fci?i -I- kp2 + kps is the same as 4,kp if 
kpi = kp4. This means that for each pinning pattern we have to carefully discuss all the ikp density operators with 
some Fermi momentum renormalizations where their wave- vectors become that of stripes. 
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We have checked through all the possible pinning patterns to see if stripes and pairing can coexist. This seems 
difficult to do since there is in principle an infinite number of pinning patterns. The fact that the pinning basis is 
determined by the orthogonal matrices Rp and Ra^ greatly reduces the difficulty. We will show that whether two 
operators can both decay with a power-law for a given pinning pattern only depends on the operator forms written 
in terms of band basis and what is the gapless mode. 

First of all, we have to know the conditions so that two (or more) operators can (or can't) be power law decaying 
at the same time. To find these conditions, consider two arbitrary operators expressed in terms of the band basis and 
let's rewrite them in the following way: 

Oa ^ ^^{^pA-Sp + '"pA-<Pp+il-aA-8a+'i'aA-<Pa) (2 22) 

Ojj ^ ^i[y-pB-Sp+VpB-$p+UaB-Sa+V„B-$a) (2 23) 

We can represent any vertex operator Oa by four coefhcient vectors UpA, VpA, u^a and Va-A- For example, if Oa ^ 
ei^{<t> i+ei) then VpA = (%/2^, 0, 0, 0), u„a = (v^, 0, 0, 0) and both UpA and v^a are (0, 0, 0, 0). 

Eq. (|2.22p and (|2.23p are written in terms of the operators in the band basis, which is not necessarily the basis 
in which boson fields are pinned. At this stage, we may not know what the new basis should be but we know at 
least that the fields in the transformed basis have to be orthogonal to each other. Then Eq. (|2.22p and (|2.23[) can be 
rewritten as 

Oa e'("p-*'^''p+"p-4''^p+'"^'i'^"+"^-^''^^), (2.24) 
Ob e*("ps-^"'p+<B-'?p+<i5-«^+»tB-'?:,)^ (2.25) 

where 6'^ — Rp9p, 9'^ — Ra^a and similarly for 0^ — Rp4>p and 0^ = Ra4'a-- Here Rp and R^r are two orthogonal 4 by 
4 matrices that transform the boson fields from band basis to a new one. The new corresponding coefhcient vectors 
are = UpARj , u'^^a ~ u^aR^^j ■■■ here _Rp and i?J are the transpose of Rp and R^. 

Now assume the bosons get pinned in this new basis. We know only one of 9[ and can be pinned and the presence 
of the dual of a pinned boson will result in exponential decay. Consider a simple situation in which all bosons are 
pinned, i.e. a COSO phase and the pinned bosons are ((/)ip, 02p> ^3p> 04p) in the charge channel and (^i^., (p^^i 4>3cr^ ^ia-) 
in the spin channel, even though we don't know explicitly the transformations to the new basis. 

This set of pinned bosons enforces some constraints on the coefficient vectors UpA/B^ "^pA/Bi '^aA/B and v^a/b so 
that Oa and Ob don't decay exponentially. The constraints are simply that the coefficients for the dual of each 
pinned boson must be zero. In this case, for X = A or B: 

^'px = (0,-,-,0), 
v'^x = (-,0,0,-), 
u',x = (-,0,0,-), 
v'^x = (0,-,-,0), 

where "— " means no constraint. These constraints actually imply the following equations: 

(2.26) 
(2.27) 
(2.28) 
(2.29) 

One can easily see that different sets of pinned bosons will imply the same equations for inner products. So the 
actually pinning patterns of boson fields are irrelevant for the constraints here and the crucial point is that all the 
fields are pinned. 

Eq. (|2.26p - (|2.29p are for the coefficient vectors in the new primed basis. Therefore, we don't really know their com- 
ponents. However, we know inner products are invariant under any orthogonal transformation. Thus, the coefficient 
vectors in the band basis should satisfy the same Eq. (|2.26p - p.29p but without the primes. So for any two operators 
Oa and Ob written in the band basis, we know the necessary condition on their coefficient vectors for Oa and Ob 
not to be exponentially decaying. 

However, things are slightly different for a CISO phase, which is the case we are really interested in. Now the 
inner products u'pj^ ■ v'^g and v'pj^ ■ u'^g in the charge channel are not necessary zero since the overlap is allowed in 

the subspace of gapless boson fields. For example none of Oa ~ e^'^'^p and Ob ^ e^^ip are exponentially decaying if 





■ v'pB 


= ^'aA 


■^aB 


= 0, 


v'pA ■ 


u'pB 


= v'aA ■ 


u'aB 


= 0, 


u'pA 


■ v'pA 


= U'^A 




= 0, 


u'pB ■ 


■ v'pB 


= u'aB 


■V'.B 


^ 0. 
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4>[p) is gapless even though • u'^g = 1 in this example. In principle, we even don't know what's the value for 
the nonzero inner products if the gapless field is arbitrary. However, here we first discuss the case in which gapless 
mode is the total charge bosons (0ip,$ip). This means, choosing (f>[p — (f>ip, that the first row of the orthogonal 
matrix Rp can be chosen to be (f/2, 1/2, 1/2, 1/2). Roughly speaking, for each charge boson field in the band basis. 



we have 



Therefore, we find that for all the pairing A and Akp density operators n^k^, the coefficients of Qip and <f>ip are fixed: 

A ~ gK'^'p&-^'p+'"'pA-$'p+K&-s',+''''a&-<l''a) _ g*-\/f *ip+---^ (2.30) 

n^^^ ~ e'("p"'*'*p+"p"''^p+"'<'"'*'*-+"--''^-) = e"'^®i''+---. (2.31) 

As we mentioned before, the only nonzero part of inner products must come from the gapless boson (0ip,$ip). 
Although we don't know the full representation of pairing and Akp density operators in the new basis, nonetheless, 
from Eq. (|2.30p and (|2.3ip . it's sufficient to conclude that the only nonzero inner product between (z/')'s and (ti')'s 
has to be exactly 



±7r, (2.32) 



where the positive sign occurs when taking the Hermitian conjugate of one of the operators. In the following conven- 
tion, we choose ikp density operators so that — tt is taken in Eq. (|2.32p . Finally, we get the necessary conditions on 
the coefficient vectors in the band basis so that a CISO phase can have both pairing and stripe correlations: 

VpA ■ Upn = -TT, (2.33) 

UpA ■ Vpn = 0, (2.34) 

"crA • Van = ^^crA ' "crn = 0, (2.35) 
UpA ■ VpA = Ua-A ■ VaA = 0, (2.36) 

Upji ■ Vpji — Ufjji • V(jji — 0. (2.37) 

What about the cases when the gapless field is not the total charge mode? We don't know the value for this 
inner product if we don't know what the gapless field is. However, we do know one thing; the inner products, Eq. 
(|2.33[) - (|2.37[) . between the coefficient vectors can only be nonzero due to the contribution from the gapless field. In 
other words, the inner products must be zero after the gapless mode is projected out. Mathematically, we mean 
'^pA ' ^pri = where v-^j^ ~ VpA — (wpA • g)g and = UpA — [upA ■ g)g where the gapless charge mode is g ■ (j>p {g is a 
unit vector). Similarly, we also have • v-^^ = 0, • v-^^ — and • = 0. The conditions on the spin fields 
are still the same since they have nothing to do with the gapless charge mode. The new conditions on the coefficient 
vectors after projecting out the gapless charge modes are 

^^A ■ ^fn = "pA ■ ^'pn 



- <A • v^n = 0, (2.38) 



^pA ■ "pA 



0. (2.39) 



Once the gapless charge mode is known, it's not difficult to check Eq. (|2.38p and (|2.39p for the pair and 4k p density 
operators. 

To illustrate what we mean by Eq. (|2.38p and (|2.39p . we take a two band system for example. The operator 
Oa ~ gi-V^Trief^+el^+tp^-i/)^) j^g^g ^jjg coefficient vectors VpA ~ \/27r(l, 1) and UpA = •\/27r(f , —1)- It's easy to see that this 
operator satisfies the condition VpA ■ UpA = 0. If (6*^^, 01^) is the gapless charge mode, then g ~ (1, l)/y/2. Only 9 
fields contain the gapless mode. The gapless mode won't appear in the cj) fields since in this case it's orthogonal to the 
(j) field, (f>i — (j>2- Then there is no contribution to the inner product from the gapless mode and we have • = 0. 
This operator could have power-law decaying correlations. However, if {9i, 4>i) is the gapless field, we have to exclude 
the part of inner product due to the gapless field since the overlap between 9 and (f> fields is allowed for the gapless 
mode. Now g = (1, 0), we find the inner product v'^a ' ^^a ~ ~27r 7^ 0. Then the correlation function of this operator 
will decay exponentially since ^2 and ^2 can't be pinned at the same time. 

It should be clear now that once we write down the bosonized expressions for pair and Akp density operators in the 
band basis, we can tell whether stripes and pairing can coexist or not by checking the inner products of coefficient 
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vectors. First consider pairing. By "pairing" we mean the existence of any operator of charge 2 whose correlation 
functions exhibit power-law decay. Any pair operator that only contains right or left fermions, such as ipRai''pRbii can 
never exhibit power law decay since from Eq. (|2.10|1 we know (pRaa = i4>aa + Qaa)l'^ and this will result in WpA • WpA 
and UctA • "SctA 7^ 0. (Eq. (|2.36p is violated.) Furthermore, any pair operator with non-zero total z-component of spin 
will contain a factor with the exponential of and hence exhibit exponential decay. 

So, there are only two types of charge 2 operators, containing only two fcrmion fields, which are candidates for 
power-law decay. Ignoring Klein factors, these are: 

i>Ra^^La^ ^ e»V^(*.P±e..) (2.40) 

^Wl.. - ) {a^h) (2.41) 

Here a and h are band indices and a =\ or |, a = —a. The -I- or — sign occurs for a =t or | respectively. Even 
though Eq. ()2.4ip carries non-zero momentum for kp^ 7^ fci?^, there is no reason to exclude it. 

We now consider the ^kp density operators, corresponding to stripes. The following two operators are the most 
general 4fci? density operators: 

- e-^^^[(^^^+^^^)+(^^r+^!:r)±(«-+«^^)±(^?r+^r)i, (2.42) 

A..^L,.A^-^^u^ - e-v^[(^^+«^^)+(^r+^- )±('^^-^^^)±(*^r-^^r)l, (2.43) 

Here «, j, k and I are arbitrary band indices and the -f or — sign is for a =| or \ respectively. Whether i, k and I 
are all different or not actually doesn't change the conclusion. When all the band indices are different, the operators 
correspond to the oscillation wave-vector of 4A:p = kpx -\- kp2 + ^fs + which is 27rn, the wave-vector of stripes. 
However, if some special renormalizations of Fermi momentum, such as kp\ = kp-j, and kpi = kp4 occur, then the 
wave vector 2{kpi + kp2) also corresponds to 27rn. Therefore, we have to consider all Akp density operators with the 
proper Fermi momentum renormalizations such that the operators correspond to stripes. That's the reason why we 
consider arbitrary rather than only all different band indices in Eq. (j2.42p and p.43p . 

In practice, we start with a pinning pattern where only the gapless charge mode is determined. Then we check 
if the pair operators Eq. (|2.40p . (|2.4ip and Akp density operators Eq. (|2.42p . (|2.43p can satisfy the inner product 
conditions (Eq. ([233]) - ((237)) if Oin.^y,) is gapless otherwise the modified Eg. ((08)) - ((09)) ). In this procedure, 
many conditions on Fermi momenta will be involved. They are the conditions to make an interaction containing a 9'^^ 
appear in the Hamiltonian (if 6'^^ is not orthogonal to 8ip, then the gapless mode is not (Oip, ^ip)), and the conditions 
to make a Akp density operator corresponding to stripes. Recall that the condition Akp = 2?™ is always satisfied. It's 
also important to keep track on the consistency of all the Fermi momentum conditions. We have carefully checked 
all possible cases and concluded that stripes and pairing can not coexist in any ClSO phases. The whole proof is 
somewhat lengthy but rather straightforward (36| . Some more details are given in Appendix A. 

Now we want to ask further if any operator of non-zero charge can have a power-law decaying correlation function 
in a stripe phase. Consider the case when (Oip, $ip) is gapless. Inspired by the real space picture of stripes and the 
finite size spectrum analysis in the next section, we find that there is always some charge-four bipairing operator, 
which does so. The most general Sz ~ non-chiral bipairing operators are: 

^Rs.^Ltr^^PRu^i',-, - e^^[(^^^+<iH(C.-+«S;)±(0:.^±0-)±('^r±«„V)]. (2.44) 

Here s, t, u and v are arbitrary band indices. Although Eq. (|2.44p are the possible bipairing operators in the most 
general sense, only when s = t and u = v or s = v and t = u docs Eq. ()2.44p carry zero momentum and have no real 
space modulation in the correlation functions. On the other hand, we have no reason to exclude the possibility that 
Eq. (|2.44p docs decay with a power law with an oscillating factor at this stage. 

Following the previous discussion, now we will prove that any ClSO phase with Akp density oscillations, also has 
bipairing correlation. The conditions for bipairing and Akp density operators to coexist are almost the same as Eq. 
p.33p - p.37p but now with Vp^hi 'Up^n = — 27r, where we use a subscript for bipairing operators. Now for simplicity, 
let's focus on the following two types of bipairing operators carrying zero momentum: 



^Rsc^Lsc^Rmi^Lm - e»v^(<*+±<r)^ (2.45) 
^Rsc^LsTi'^Rtii^Ltp - e'^''^'st±e:t)_ (2.46) 

We find that Eq. (|2.42p . (|2.43p and Eq. (|2.45p can coexist if we choose {s, t} = or {/c, /}. If we have a = (3 in 

Eq. (I23S1), i-e. 6*^+ is present, then Eq. ^^T^ . ^(I^ and Eq. ^(XWf can coexist with the choice that {s,t) = {i, j} 
or {fc, I}. If we have a — (3 m Eq. (|2.46p . then 91^^ is present. Eq. (|2.42p and (|2.46p can coexist with the choice that 
{s, t) = {i, k) or {j, I). Eq. ((^15)) and ^(XWf can coexist if {s, t} = {i, 1} or {j, k). 
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C. Completely Empty or Filled Bands 

Another possible type of stripe phase has one or more of the bands completely empty (or completely filled). We 
now show that a standard treatment of these phases does not lead to any with coexisting pairing and stripes. Suppose 
that two bands are completely empty. Without loss of generality, we may choose then to be bands 3 and 4 so that all 
the electrons go into bands 1 and 2. It then follows from Eq. (|1.7[) that 

fcFi + fcF2 = 27rn = 27r(l - 5). (2.47) 

So stripes, i.e. Fricdcl oscillations at wave-vector 2ttti, corresponds to 2kp oscillations in the effective 2-band model. 
There are two cases to consider: 

A^^m. - e^^(«?±'^n (2.48) 
V'k^fi,. - e'v^(^S'-^^7±^?^T0r). (2.49) 

For the above situation, we have {i,j} — {1, 2} in Eq. (|2.48p and (|2.49p . Following the similar method, we can prove 

the incompatibility of pairing and 2kp oscillations. If instead bands 3 and 4 are completely filled, Eq. (|1.7|) now 
implies: 

kpi + kF2 + 27r = 27rn, (2.50) 

but this is equivalent to the case of the two bands being empty since exp(27rmj:) ~ exp[27ri(n— l)x] for any lattice site, 
X. If only one band (4) is empty (or filled) then stripes would correspond to oscillations at wave- vector kpi + kp2 + kF3 
which can never occur since any operator which occurs in the continuum representation of the density operator must 
contain an even number of fermion fields. Similarly stripes could not occur in a phase with 3 empty (or filled) bands 
since this would require oscillations at kpi. 



D. Summary and Generalization 



Therefore, for the CISO phase with stripes, there is no pairing but bipairing occurs. In Sec. V we will discuss 
the case of small t_L,2, Vj_,2 and the analogy with a 2-leg bosonic ladder. The bipairing correlation in the four-band 
fermion system may correspond to the boson pair superfluid (BPSF) phase in Ref. (27| . 

Now we would like to generalize this argument a little bit further. What about any other generalized charge two 
operators such as 

^Ha.i^L,.i^l^^H,. - e'^^[('^- +^^) + (^«ir-0-)±K^+^-)±(0:^-0?r)l. (2.51) 

Although the correlation amplitude should be smaller compared to the usual pairing operators, still, is it possible that 
such generalized pairing operators after renormalized by some density-like (charge neutral) operators can coexist with 
stripes? There will be a lot more such charge two operators since Eq. (|2.40p . (|2.4ip and other unconventional pair 
operators can be combined with any charge neutral operator, as long as in the end we have a non-chiral, spin zero 
and charge two operator. So far we have checked up to charge two operators composed of six fermions and concluded 
none of them can coexist with stripes as expected. But what about charge two operators composed of more fermions? 
This endless question may require another approach for its resolution. Instead, we will study the finite size spectrum 
of a CISO phase, from which a connection between charge operator and the lowest density oscillation is established. 



III. FINITE SIZE SPECTRUM 



In this section we will establish a general connection between charge and density operators through the consistency 
for the finite size spectrum of a CISO phase. Assume that the total charge field is the gapless charge mode. The low 
energy effective Hamiltonian in a CISO phase is simply that of a free boson. 



H 



Ip 



(3.1) 



In the rest of this subsection we only discuss the boson <I>ip (and its conjugate boson Qip) so, for convenience, in this 
section we drop the superscript p and the subscript 1. Here vdxO = Kdt^ = 11, where 11 is the canonical momentum 
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variable conjugate to It is natural to regard $ as a periodic variable. Consider first a pairing phase where a charge 
two operator: 



A{x) - eW2^^ (3.2) 

has power law decay. Only keeping operators in the low energy Hilbert space which exhibit power-law decay . we 
expect that all such operators will have even charge, involving only the exponential of integer multiples of It 
is then natural to assume that we should make the periodic identification: 

$ ^ $ + 2V2^. (3.3) 

That is to say, we regard •y7r/2$ as an angular variable. We now wish to argue that consistent quantization of the 
free boson requires that O also be regarded as a periodic variable with: 

e e + y|. (3.4) 

As we shall see, this, in turn implies that the minimum Friedel oscillation wavevector is Skp. Alternatively, in a 
bipairing phase, the lowest dimension charge operators is exp[i-\/27r$J] and it is now natural to identify 

$ ^ $ + \/2^. (3.5) 

We then will argue that consistent quantization requires: 

e ^ e + \/2^, (3.6) 

which we will show implies that the minimum Friedel oscillation wavevector is Akp. This approach confirms the 
conclusions arrive at by more pedestrian means in the previous sub-section. As a biproduct of this discussion, we will 
derive the finite size spectrum, with both periodic and open BCs, in both pairing and bipairing (stripes) phases. 

First consider a CISO pairing phase. We place the system in a box of size L with periodic boundary conditions. 
The mode expansion for $(t,x) takes the form: 



' k=i » ^ 



(3.7) 



Here m and p are arbitrary integers; a^k and a^k are bosonic annihilation operators for right and left movers. The 
normalization of the mx/ L term in the expansion is determined by the periodic BCs and the identification Eq. (|3.3p . 
i.e. ^{L) ~ $(0)-|-2V27rm is equivalent to PEC using Eq. (|3.3p . The (very important) normalization of the vt/ L term 
requires more explanation. We may think of this term as being proportional to a zero mode conjugate momentum 
operator Ho, which is canonically conjugate to <I>o' 

[$o,no]=i, (3.8) 

$(t,x) = $o + ^^ + ... (3.9) 
KL 

Ho is the zero momentum Fourier mode of the conjugate momentum field n(x): 

Ho = / n(x). (3.10) 



JO 

The T\f)vt/{KL) term in the mode expansion for $ is necessary in order that the canonical commutation relations are 
obeyed: 

oo 

[$(x),9t$(y)] = -^ ^ e*(2-'=/^)(-.)^^5p(.^_y), (3.11) 

k— — oo 
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where 5p(x — y) is the periodic Dirac 5-function. The k = Q term in this Fourier expansion of the (^-function comes 
from the commutator in Eq. p.Sp . Comparing the mode expansion in Eq. (|3.7p to p.9p . we see that the eigenvalues 
of the canonical momentum operator, Ho are: 



(3.12) 



no = \ ji^p, 



for integer p. That these are the correctly normalized eigenvalues follows from the fact that the wave-functionals 
contain factors of the form: 



*(*o) 



ino*o 



(3.13) 



These wave- functions are single- valued under the identification of Eq. (|3.3p . The eigenvalues of Ho, i.e. the normal- 
ization of the vt/ L term in the mode expansion of Eq. (|3.7p determines the charge quantum numbers of all low- lying 
states in the spectrum with periodic boundary conditions. This follows from observing that the total electron number 
operator is: 



2K 

V 



\ - I dxdt^ 



TT Jo 



(3.14) 



This in turn can be checked by confirming that: 



where A is the charge-2 operator in Eq. 
states in the low energy spectrum are: 



[7Ve,A]=2A, (3.15) 
The mode expansion of Eq. p.7p then implies that the charges of all 



Ne = 2p, 



(3.16) 



even integers. This is an example of the general one-to-one correspondence between operators and states in the finite 
size spectrum with PBC in a conformal field theory. In a phase in which all operators have even charge, all states in 
the spectrum also have even charge. 

We can go further and deduce the Friedel oscillation wavevector from the normalization of the vt/L term in Eq. 
p.7|) . This can be done by using Kdt^ = vdxQ, dtQ = vKdx^, to deduce the mode expansion for the field Q{t,x): 



Q{t,x) = eo + J-Pj- 



2\/2TTKm- 



vt 



E 

fc=i 



K 
47rfc 



From this mode expansion we see that 9 is periodically identified as in Eq. (|3.4p . Insert Eq. p.7p and (|3.17p into the 
Hamiltonian (|3.ip . we obtain the finite size spectrum of a pairing phase with PBC: 



E-2fip 



2ttv 



2Km^ 



8K 



k=l 



(3.18) 



where n^^ and njik are the occupation numbers for the left and right moving states of momentum -izlirk/ L. 

Now consider the Friedel oscillations. Oscillations at wave- vector of the form 2nkp (actually a sum of any 2n of the 
kpi's) can only occur if some operator of the form {ippipL)" has power law decay. Upon bosonizing, all such operators 
are expressed as exp (^in^/^T/2Ql^ multiplied by an exponential involving only pinned boson fields. However, not all 

such operators can occur in the low energy spectrum since they must respect the periodic identification in Eq. (|3.4p . 
The lowest dimension operator allowed by this identification has n = 4 corresponding to Skp oscillations. Again we 
are effectively using the relationship between the finite size spectrum and the operator content. The ?i = 4 operator 
corresponds to the p = I state in the mode expansion of Eq. (|3.17p . 

Now consider a bipairing phase where there are no charge 2 operators in the low energy spectrum, the lowest charge 
being 4, corresponding to the operator exp[i\/27r$], leading to the periodicity condition on <I> in Eq. (|3.5p . The mode 
expansion for $ is therefore altered to: 



Ht,x) 



27rm- 



/27r 



p vt 



E 

k=l 



1 

AirkK 
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The coefficient of vt/L gets multiplied by a factor of 2 since the wave-functional exp[ino<I>o] must now be invariant 
under the shift of Eq. p.Sp . requiring the conjugate momentum, Ho to have eigenvalues •\/27rp. Correspondingly the 
mode expansion for becomes: 



A;=l 



implying the periodic identification of Eq. (j3.6p . Now the lowest dimension Friedel oscillation operator with power 
law decay is exp(i'\/27r8]'), which is a Akp operator. The finite size spectrum of a bipairing phase with PBC is: 



E ~ 



2TTV 

~17 



9 OO 



fe=l 



(3.21) 



These arguments show, based only on plausible assumptions about regarding the fields <i>j and 8j as periodic 
variables, that CI SO phases with pairing have 8/cf oscillations (and hence no stripes) but phases with bipairing have 
oscillations, corresponding to stripes. 
With OBC, the boundary conditions: 

6(0) = constant, Q{V) = constant, (3.22) 

are applied. This sets the quantum number m = and a_Rfc = ULk in the mode expansion of Eq. (j3.17p and (j3.20p . 
Setting Q = 2p for the pairing phase and Q = 4p for the biparing phase, the finite size spectrum implied by these 
mode expansions and the Haniiltonian of Eq. (|3.ip is: 



WK 



k=l 



(3.23) 



where the charge, Q (measured from a reference point like half-filling) is restricted to all even integers in a pairing 
phase but is restricted to integer multiples of 4 in a bipairing phase. For even Q, in a bipairing phase, there is a gap 
A^, to states with Q/2 odd, so we may write, for any even Q in a bipairing phase: 

fl- (-1)^/21 I 
E-fiQ = AE^- ^ + (3.24) 

The parameters v, K, fi and all vary with density. Nonetheless, this zigzag pattern of energies for even Q should 
allow unambiguous detection of a biparing phase for large enough L. 

IV. RG FOR DOPED 4-LEG LADDER 

The combination of weak-coupling RG and bosonization is one standard tool to study the phase diagram of TV-leg 
systems [21,22. 26,.37| . The results for doped 4- leg ladder are mostly within the context of the Hubbard model 2l|[2^ . 



Here we will show that the stripe phase can be found in the special solution of RG equations. This phase doesn't 
have pairing but bipairing, which is consistent with our bosonization argument. 

The first step is to determine the relevant couplings according to the RG flow since they will control which boson 
fields will get "pinned" and therefore will allow us to map out the phase diagram in terms of bare interactions and 
doping. However, to analyze the RG flow is a tricky task for there are 32 coupled nonlinear differential equations. 
It seems that the RG ultimately flows onto a special set of solutions, corresponding to some direction in the multi- 
dimensional coupling constant space. These special solutions are called "fixed rays" and different rays usually indicate 
different phases plj . Some fixed-ray solutions may correspond to phases with higher symmetry than the original 
Hamiltonian. Two-leg ladders at half-filling provide one example of such symmetry enhancement in the low energy 
limit [38| . Later on, some subleading corrections were found which make the RG flow deviate from the flxed ray. 
However, these subleading terms don't grow fast enough to really spoil the fixed ray in the undoped case but give 
some anomalous corrections to the gap functions, vanishing in the weak coupling limit 2^|3^. Things become 



dramatically different in the doped systems. Now these subleading terms are relevant perturbations for the fragile 
gapless modes. They will generate gaps although these may be much smaller compared to those driven by the "fixed 
ray". For example, the weak coupling RG phase diagram for doped two-leg ladders is modified after taking these 
terms into consideration [2^. Recently, a hidden potential structure of RG equations in ladder systems was discovered. 
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Everything can be understood better within the framework of this "RG potential" which allows the RG flow to be 
viewed as the trajectory of a particle finding a minimum in the coupling constant space (2^l25j . Then the fixed ray 
is just like a "valley/ridge" in the "mountains" of the RG potential. The topography near the vicinity of such a 
"valley/ridge" will determine the stability of the fixed ray and give the exponents governing the subleading terms. 



A. RG Potential and Its Implications 



The method we discuss is very general but we mainly focus on iV = 4. As the conventional starting point, we first 
diagonalize the hopping terms in Eq. (|2.ip and obtain 4 bands. Next wc linearize each band around different Fermi 
points in the low energy limit, and introduce SU{2) scalar and vector current operators 



J^/^ = ^ V'L/flja , Jf/^ = \ i^l)R, i'L/RjP ■ (4.1) 

(Note the unconventional factor of 1/2 in the scalar operators, introduced for later convenience.) We can rewrite the 
interactions in Eq. (|2.2|) in terms of the current operators in Eq. (|4.ip : 



where fij and denote the forward and Cooper scattering amplitudes, respectively, between bands i and j. Many 
repeated indices appear in this section, such as i and j in Eq. (|4.2p and they are always implicitly summed over. In 
Eq. (|4.2p . we only keep the Lorcntz invariant interactions involving the product of a left current and a right current. 
The LL and RR terms don't contribute to the RG equations at second order and are expected to only shift the 
velocities of the various modes. Note that cu and fa describe the same vertex so we set fa = 0. Also, symmetries 
imply Cij = Cji and = fji [2l|. That's how we get 32 different couplings in doped 4-leg ladders. Then one can 
derive RG equations by the operator product expansions of these SU(2) scalar and vector current operators. 

Provided that the RG equations arc known [2l|, how to analyze the RG flow is still non-trivial. Since all the 
interactions in Eq. (j4.2p are marginal at first glance, if one numbers all the couplings fij and and rename them 
as gi {i from 1 to 32), within the one-loop calculations, the coupled non-linear RG equations can be written in the 
concise form: 

^ = Mtmk, (4.3) 

where gi is some coupling and the coefficient matrices M^^ = M^*'-' are symmetric in indices j and k by construction. [I 
is the logarithm of the ratio of a characteristic length scale to the lattice scale, I = ln(i/a).] Recently, an unexpected 
potential structure of Eq. ()4.3p was proven (23l25l |. After a proper rescaling to new couplings gi = aigi, where ai are 
constants, Eq. (|4.3p can be reduced to 

where Mil'' is totally symmetric in indices i,j and k and V{g) is the so called RG potential [2^. The scaling constants 
ai and the explicit RG potential form can be found in Appendix C. It now provides a geometric picture for the 
RG flows of Eq. (|4.4p . which can be regarded as the trajectory of an overdamped particle searching for a potential 
minimum in the multi-dimensional coupling space. Thus, the ultimate fate of the flow would either rest on the fixed 
points or flow along some directions as the "valleys/ridges" of the potential proflle but there is only a trivial flxed 
point (all gi = 0) within one-loop order. 

Precisely, these directions are special sets of analytic solutions of Eq. (|4.4p : 

5.(0 = 7^, (4.5) 

I'd ~ t 

if the constants Gi satisfy the algebraic constraint, 

G, = A/f GjGfc. (4.6) 

Eq. (|4.5p is only valid ioi I < Id = In^/a where ^ is a characteristic length scale where the coupling constants become 
large. These special analytic solutions are referred as "flxed rays" because they grow under RG with the fixed ratios. 
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Sometimes, the specific ratios of the fixed ray reflect extra symmetry in the Hamiltonian and the flxed ray is called a 
"symmetric ray" [38j . 

In general, it's very unlikely that the bare values of gi{0) are proportional to the constants G;. Then we have to 
check whether these fixed rays are stable against deviations [25|. As long as the deviations grow slower than Eq. 
(|4.5[) . then the fixed ray is stable. Within the RG potential picture, the stability of each fixed ray is determined by 
the local topography along the direction. 

In the vicinity of the fixed ray, if there are some small deviations away from it 

g^{l) = + Ag^il), (4.7) 

where Agi{l) <C gi{l)- The equations which describe the deviations Agi are 

= A,, , ,4.8) 

where By — 2A/j^Gfc. Since Mjf. is totally symmetric in i,j and k, the matrix Bij is symmetric in i and j. Bij can 
be diagonalized by an orthogonal matrix Onm so that Eq. (|4.8p will decouple into independent equations. 



|(^.„) -J^^Sg.^ (4.9) 

where Sgn = Oni Agi, are the couplings after the linear transformation and A„ are the eigenvalues of the matrix Bij. 
If the initial bare couplings 6gn{0) ^ Gn/ld, then the solutions of Eq. (|4.9|) are 

Sg.(l)^Sg,m{j^^)'^ (4.10) 

where G'„ ~ 0{{U/ty~^"') is generally non-universal depending on the initial couplings. Therefore, the appropriate 
ansatz for the RG flows should be 

- ■ ^' +•••, (4.12) 



Id -I {Id- 0^'°°" 

where Af^'^ in Eq. (|4.12[) is the largest one among A„'s with nonzero coefficients OmGj^ = G" and the divergent 
behavior is dominated by A™'*'^. Although Eq. ()4.12p is derived from the stability analysis near the fixed ray, it seems 
to be the general behavior of the RG flow but the values of some A™'^'^ may vary from the eigenvalues of Bij when 
away from the flxed ray. In fact, such power-law divergent solutions were suggested before |2lll37| but the analysis 
only focused on the most relevant terms with exponent one, i.e. Gi ^ terms. Note that Eq. (|4.12p is still not the 
exact solutions of RG equations but it captures the divergent part correctly and is enough to determine the phase 
diagram. 

For A™'''^ < 0, these deviations are irrelevant whereas if A™ '^'^ > 1, the deviations grow faster than the fixed ray, so 
that the phase associated with the fixed ray is fragile. For < X™^^ < 1, the deviations actually grow although not 
strongly enough to spoil the asymptotic fixed ray, as illustrated in Fig. ^ . Nevertheless, it doesn't mean they won't 
affect anything since they are also relevant couplings in the conventional classification but just less relevant than the 
fixed ray ones. 

The effects of these subleading divergent terms on the RG flow of a particular coupling gi are dramatically different 
depending on whether G, = or not. We can separate the effective Hamiltonian at the cutoff length scale into the 
most relevant fixed ray part and the subleading deviations as a perturbation. For Gi ^ 0, gi is relevant and it will 
lead to a gap. The subleading perturbations will not destroy the original ground state but only modify the gap 
function by giving rise to anomalous scaling [2^ . On the other hand, if Gi = 0, the subleading perturbations become 
important if < A™'^'^ and will generate a small but non-zero gap out of the initial gapless modes. Note that only 
when initial deviations away from the fixed ray are small, are Xf^^^ universal and can be obtained from the eigenvalues 
of Bij = 2Mjj,Gfe. However, from numerical solutions of the RG equations we find that the couplings always diverge 
with power law behavior like Eq. (|4.12p though Xf^'^^ is not the same as the eigenvalues of Bij . We can extract Xf^'^^ 
directly from the numerical solution of the RG equations for those terms with Gi = and thus determine the phase 
diagram from these relevant interactions. Surprisingly, as we will see in the following sub-sections, the fixed ray and 
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l<0 0<A,<1 I >l 



FIG. 1: The topography of RG flows near the fixed ray with A<0, 0<A<1 and A > 1. It's clear that the deviation is 
irrelevant for A < and relevant for A > 1. The analysis of RG flow is more subtle for < A < 1. In this case, although the 
deviation from fixed ray is growing, the fixed ratios still remains. Therefore, RG still flows onto the fixed ray but the phase is 
not only determined by the fixed ray couplings. 

subleading terms with universal Xf^'^^ from Bij already are sufRcient to determine the pinned bosons uniquely within 
bosonization. 

To recap, even in the weak-coupling RG analysis, there are different energy scales. The fixed ray only represents 
the most relevant couplings. The subleading couplings should be treated as perturbations to the effective Hamlitonian 
corresponding to the fixed ray and they are relevant enough to drive the effective Hamiltonian into a phase in which 
the fragile gapless modes become gapped but these gaps are small compared to those driven by the fixed ray couplings. 



B. The Stripe Phase 



Now we know the ultimate fate of weak coupling RG must be a fixed ray due to the existence of the RG potential 
l^H^. The fixed ray indicates a direction in which the interactions will be renormalized in the strong coupling region. 



If one can survey all the fixed ray solutions in Eq. (j4.6p and the corresponding subleading terms determined by the 
topography, then in principle, all the phases in the weak coupling RG are obtained. Following this idea, here we try 
to find the fixed ray whose corresponding phase gives the stripe density oscillations. It turns out that following fixed 
ray will do so: 



2 2V3 2V3 

2 2V3 ^ 2x/3 16{ld-l) 

Eq. (|4.13p will be the solution of Eq. (|4.6p if I'l = V4 and V2 = v^. So if RG really flows onto this fixed ray from some 
initial set of bare couplings, the interpretation is that the fermi velocities get renormalized in the corresponding phase 
[sst . The upper and lower line in Eq. (|4.13p correspond to the CDW fixed ray on effective 2-leg systems composed of 
band pairs (1,4) and (2,3), respectively. In principle, the fixed ray as the permutation of band indices in Eq. (|4.13p 
also exists. The reason why we favor Eq. (|4.13p is motivated by the fixed ray Eq. (|4.19p . found in the Hubbard model 
as we will see in the later sections. 

Now we know that the fixed ray solution isn't the whole story for the RG flow. The phase should be determined 
by all the relevant interactions, including the subleading divergent ones. As long as Gi is given by the Eq. (|4.13p 
and with the known RG matrix M/'^, then the largest divergent exponent X™'^'^ can be deduced analytically from the 
eigenvalues and eigenvectors of the matrix Bij = 2Mjj.G'fc in the vicinity of the fixed ray. There are two 5/8, four 
1/2, four 1/8, four —1/2, two —3/8, and Os for the eigenvalues A„ of Bjj. The couplings are divergent for A„ > 
even though they are small compared to the fixed ray in the critical region. We should take more care about the 
terms with eigenvalues 5/8, 1/2 and 1/8. To see what's their influence, we have to know the direction corresponding 
to 6gm Eq. (j4.10p . in the multi-dimensional space expanded in the coupling basis Agi. The eigenvectors of Bij 
give this information. The subleading terms dgn corresponding to two A„ ~ 5/8 eigenvectors are in the directions 
having non-zero projection on c^2, c^g, C24, C34, Ci2-, cjg, 0^4, and C34. The eigenvectors of those corresponding to four 
A„ = 1/2 have components on c^i, c^2' C33, 0^4, ^4, 4^, Cn, 0^2^ C33, 0^4, c^^, cjg, f^^, f^^, J^^ and /2''3. The terms with 



18 



A„ = 1/8 have components on c^2 1 '^13 ' "^24 > '-34 ' ^12 j c^3 j C24 , and C34. Provided with this information, we know the 
largest divergent exponent for each couphng, that is, X'"^'^^ for the non-fixed ray coupHngs. Table ([I]) summarizes A„ 
and the projections of corresponding Sgn in terms of couphng basis gi . 





nonzero component 


^ /ft 

0/0 


„P „P „P „P ^O" r-" 
''121 '-13' '-24' '-34 1 '-12' '-34 


^ /s 

0/0 


nP nP r-P n'^ n'^ 
^-12 ' '-IS ' '-24 ' '-34 ' '-13 ' '-24 


1 /9 


rP r-P rP fP f 
'-11' '-44' '-14' '-11' '-44' ./l4' ./14 


1/2 


^P r-P r-P n" r" fP f 
'-22 ' '-33 ' '-23 ' '-22 ' '^33 ' J 23 ' J 23 


1/2 


r-P r" r" r" 
'-14' "-ll' <-44' '-14 


1/2 


(^23 ' "^22 ' <333 ' C23 


1/8 


^P „P 
'-12' '-34 


1/8 


„P „P 
'-13' '-24 


1/8 


^12 ' C34 


1/8 


(T cr 
Cl3' '-^24 



TABLE I: This table summarizes the topography in the vicinity of the fixed ray Ea. (|4.13p . It shows the eigenvalues \n > of 
the matrix Bij and their corresponding eigen-direction in terms of the RG couphngs. 



Now we would like to check numerically whether RG will really flow onto this fixed ray Eq. (|4.13p . It's very 
illuminating to plot log[|gi(?)|] v.s. log[(/d — I)] from the numerical solution of the RG equations, where the absolute 
value makes sure there won't be problems for those with Gi < 0. In the scaling region, if Eq. (|4.12p is correct, then 
we should see a straight line for each coupling gi{l), whose slope indicates the exponent controlling the divergence. 
The slopes will be negative one for the fixed ray, —X^'^^ for the subleading terms and zero for irrelevant terms. 

If we choose the initial bare couplings with the ratios in Eq. (|4.13p and with Fermi velocities vi = ^4 and V2 = 1^3, 
we do find all the couplings grow with the fixed ratios under RG flow toward the fixed ray Eq. (|4.13|1 . In order to see 
the subleading terms, wc add some small deviations to the initial bare couplings and Fermi velocity. The log-log plot 
of each coupling agrees very well with the predicted slopes in Table (|l|. A few selected examples arc shown in Fig. ([2]). 
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FIG. 2: This is the log[|(ji(?)|] v.s. \og[{ld — I)] plot for several typical couplings of stripe fixed ray Eq. (|4.13|) . The slopes give 
us the divergent exponent of each coupling. The solid (red) lines are the numerical solutions of the RG equations. The dashed 
lines are pure straight lines as reference with the predicted A™^": 1 (pink: C14 (a), /fs (b)), 5/8 (blue: b^^ (c), C34 (d)), 1/2 
(green: c^i (e)), and (yeUow: /fs (f)), respectively. In this case, the numerical solutions agree very well with the prediction 
for aU the couplings. 
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Although the subleading terms are also divergent for < X™^^ < 1, they are still small compared to the fixed ray 
couplings. Therefore, we treat the subleading couplings as the perturbations to the effective Hamiltonian corresponding 
to Eq. (|4.13p in the bosonization method. 

We bosonize the relevant couplings to determine the phase diagram. The full bosonized form of Eq. (|4.2p is 



Hi, 



4 « « 



cScos(V87r6',^) 



J Zlt^^y- + %/4^<" cos Vl^e^- + 2dfj cos V4^(f>^~ cos Vi^e^+ 



+ (C - cfAcos V4^0f- cos \/4^0^- - 2/^- cos \/4^0°." cos VI^0: 



-+1 
ij J 



(4.14) 



where these bosons fields are defined in Eq. (|2.14p . Since cf^ and f[j only contribute to the gradient terms, they are 

not important here. The reason we express the Hamiltonian by the tilde interactions c, / and Ea. (j2.14p is only for 
the convenience of notations. We emphasize that the basis of boson fields should be determined by the hierarchy of 
relevant interactions. In other words, the relevant interactions should not only tell us what boson fields are pinned 
but also the basis in terms of which they are pinned. 

At some intermediate length scale, the most relevant interactions, those in Eq. (|4.13p . are large but the others, 
including the subleading terms, are still small compared to them. In order to minimize Eq. (|4.14p . the most relevant 
couphngs in Eq. (|4.13p will pin the values of 't'u t 4>23 ^ ^it ^ ^23 ^ 'I'i4 ^^'^ 4'23 ■ These pinned bosons, regardless of 9 
or 0, immediately suggest a basis. In the spin channel, since there are already four mutually orthogonal combinations 
of band bosons getting pinned, it's natural to choose Ra corresponding to the pinned combinations. So the basis 
for spin fields is fixed. As for the charge channel, there are two combinations of band bosons which get pinned and 
we also know symmetry requires the gapless mode to be total charge field ($ip, 6ip) since there is no interaction 
involving boson fields not orthogonal to it here. These three fields and the orthonormal condition will uniquely fix 
the only unknown basis field. 



1 



2p 



hp - 4>2p - hp + <t>ip), 



(4.15) 



similarly for 02p if replace 
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by 9. In this case, the relevant interactions suggest the basis we should adopt is: 
1/2 \ 
-I/V2 
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Now we will rewrite the interactions Eq. (|4.14p in terms of the new basis given by Rp and Ra- So far we get six 
pinned bosons only by considering the fixed ray interactions. At lower energy scale, the subleading terms also become 
large, yet still small compared with the fixed ray. Replacing the pinned bosons by their pinned values, we get a C2S0 
effective Hamiltonian and the subleading terms will be treated as the perturbations. $ip is absent in Eq. (|4.14p and 
it will remain gapless as we expected. The question is about whether $2p will get pinned due to the perturbation 
involving it, such as c^2i ^13, C24, C34, c^2, C24 and cf4. At first glance, one may conclude the gapless bosons <I>2p 
won't get pinned since all the subleading perturbations also contain the dual of pinned spin boson 0^4" ^'^d 023~ ■ 
other words, these subleading interactions should be irrelevant and <I>2p can't be pinned. It's true for this analysis. 
But the common wisdom tells us that the gapless mode is usually fragile unless protected by some symmetry or 
incommensurability. In fact, there are always some other higher order interactions which can be generated in the 
continuum limit as long as they are allowed by symmetry. We usually don't pay attention to these higher order terms 
for they should be much smaller and less relevant than the interactions in Eq. (|4.2p . However, the scaling dimension 
of these higher order interactions can be changed due to the existence of some other interactions (39ll40| . For example, 
consider a 4*'' order term in perturbation theory: 



SH (X (c?2) 



[cos V A-K (1)^2 COsV2tt{9i^ + 92a)] cos VStT^Ict COS V 87r02cr 



(4.16) 



Using the operator product expansion, we can replace all factors involving 9i^ and 92a by a constant. The remaining 
operator contains a term: 



5H cx cos v87r( 



^ip 



'^2p. 



= cos\/47r(%/2$ 



^2p 



i?4 



^'23) 



2p, 



(4.17) 



where we replaced cj)]^^ and 023 by their expectation values in the third expression. This operator doesn't depend on 
spin fields anymore and we have an effective sine-Gordon Hamilatonian for ($2p,02p)- The cosine interaction has a 
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scaling dimension of IjK^p- If the rcnormalizcd value of the Luttinger parameter for the <I>2p boson, 7^2p > 1, then 
Eq. (j4.17p is relevant. In general, it's highly nontrivial to determine the renormalized value after integrating out 
the gapped modes. However, we can calculate the rcnomalization of due to the gradient terms of interactions: 

K^p^J-^-^, (4.18) 

y TTV + C — J 

where 

V = {Vi + V2 + V3 + V4^) / 4, 
C = C^l + + C33 + C44, 

/ = .fl2 + .fl3 ^ flA ~ /23 + /m + /34- 

According to the ratios in Eq. (|4.13p . we find that K2p > 1. Therefore, $2p should also be pinned with this sine-Gordon 
type interactions. We conclude the final phase should be CI SO and the pinned bosons are (f>2p, 014^, ^23^5 ^14^7 ^23'' 't'14 
and cf)2^ . For this pinning pattern, the correlation functions of the 4fc p density operators Eq. (|2.42p and (|2.43p with 
{hj} = {lj4} and {k,l} = {2,3}, decay with a power-law. Also, with {s,t} = {1,4} or {2,3}, the correlation 
functions of the bipairing operators Eq. (|2.45p and the term with O^j^ in Eq. (|2.46p decay with a power-law. All the 
pairing operators Eq. (|2.40p and (|2.4ip decay exponentially. This phase has stripes and bipairing correlations. The 
result is also consistent with FSS. 

The only question left is how to find the proper initial couplings so that the RG will fiow to the fixed ray. The 
initial bare couplings are determined by the interactions in the model. As long as one includes enough short ranged 
interactions in the Hamiltonian, the initial bare couplings can be tuned near the ratios in Eq. (|4.13p . The point is 
that the fixed ray should indicate some phase in the strong coupling regions. Thus, to find the proper initial bare 
couplings that RG will fiow to this fixed ray may not be the most important issue for our purpose. 

The fact that we need vi ~ V4 and V2 — in order to see the stripe phase in RG resembles the situation in the 
decoupled 2-leg ladders limit we study in section V. It seems to suggest that the stripe phase should be related to the 
renormalization of Fermi velocities from both limits we study. 



C. The Weak Coupling Repulsive Hubbard Model 

In the previous section, we found the fixed ray corresponding to the stripe phase without knowing the exact 
underlying model. Now we would like to switch gears and study the fixed rays corresponding to parameters of the 
Hubbard model in Eq. ^Ji) and (j^ . 

In the four-band region, the RG flows to the following fixed ray (2l| . 

= = 4 . 4f . . . -/f. . -j^. (4.1,) 

It seems that only the interactions between band 1 and 4 are relevant and band 2 and 3 are totally decoupled from 
the system. Also, the fixed ratios in Eq. (|4.19p are the same as those of the CISO phase in a doped 2-leg ladder if 
band 1 and 4 are regarded as an effective 2-leg system. So the final phase was referred to as CISO -I- C2S2 = C3S2. 
Notice that the RG equations, in the gi basis of Eq. (|4.4p . are invariant under the permutations of indices but the 
bare values, 5i(0), favor the phase in which the couplings involving bands 1 and 4 get large. This can be seen from 
the factors of 1/vi relating the g^'s to the g^'s in Eq. (|C9P - (jC12p . 

However, now we know that the fixed ray solution isn't the whole story for the RG flow. Bands 2 and 3 are 
never really decoupled since there are subleading coupling constants that involve these bands. Once again, we will 
plot log[|gi(Z)|] v.s. log[{ld — I)] from the numerical solution of the RG equations and we see nice straight lines in 
the scaling region. The slopes will be compared with —A™'''', which can be deduced from the eigenvalues A„ and 
eigenvectors of the matrix Bij = 2M^-^.Gk- 

With the Gi given by Eq. (|4.19p . there are two —1/2, six —1/16, two 1/2, two 15/16, a 1 and other Os for A„. 
Terms corresponding to —1/2 and —1/16 are irrelevant and not harmful to anything. We should carefully look at 
the terms with 1/2, 15/16, and 1. The effects on the phase diagram depend on whether they have components on 
the couplings with = 0. This information is given by their corresponding eigenvectors. One may think in general 
A„ = 1 means the fixed ray is unstable. This is true if the deviations Sg^ have non-zero components on the couplings 
besides the fixed ray ones. Fortunately, here the eigenvector for A„ = 1 only has two components with negative c^^ 
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and positive C44 in terms of the original coupling basis. It will shift the values of fixed ratios regarding c^^ and C44 in 
Eq. (|4.19p a little bit but it won't change the fact that those seven couplings are the most relevant ones. This only 
reflects that the fixed ray Eq. (|4.5p . is just a special set of the solutions and not the most general one. Table pl|) 
summarizes A„ and their projections in terms of coupling basis gi. 
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p p p a a a rp 
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TABLE II: This table summarizes the topography in the vicinity of the fixed ray Eq. (|4.19|) . It shows the eigenvalues A„ = 1 of 
the matrix Bij and their corresponding eigen-direction in terms of the RG couplings. 

This result can be checked by plotting log[|gi(Z)|] v.s. log[{ld — I)] for the numerical solutions of RG equations. As a 
test, we can artificially tune the ratios of initial conditions based on Eq. (|4.19|) . such that the RG flow will be really in 
the vicinity of the fixed ray. In this case, the slope of each coupling agrees perfectly with the prediction given above. 

Now using the initial conditions as shown in Appendix C, determined by physical parameters on-site interaction 
U and the doping S, we can do the same analysis. Even in the region RG flow still controlled by the flxed ray Eq. 
(|4.19p . now we find not all the slopes agree with the prediction. Some couplings with predicted Xf^^^ = 0, actually 
have non-zero slopes in the log-log plot and those slopes may vary according to the initial conditions. They are new 
subleading terms besides those given by the stability check near the fixed ray. However, we don't find the notable 
change of the slopes for those couplings with \™^^ — 5/16, 1/2 or 1, as long as the RG flow is still dominated by the 
same fixed ray. That is, the universal analytic prediction in the vicinity of the fixed ray is still correct to some extent. 

A few selected typical examples are shown in Fig. (jS]). The failure to predict all Af^'^ correctly for each couplings 
doesn't mean we can't determine the phase diagram. The point is that we should treat all the divergent terms with 
the exponent < Xf^'^^ < 1 as the perturbations to the effective Hamiltonian corresponding to the fixed ray in the 
bosonization scheme. 

In order to minimize Eq. (|4.14[) , the most relevant couplings in Eq. (|4.19[) will pin the values of (jy^^ , 9^^ , and 
. The fixed ray interactions and the symmetry imply we should choose $2' 4>i4j (^14 j ^i4~ ^^r the new basis 
of bosons. Unlike the previous case of stripe fixed ray, here there are still two undetermined fields in charge and spin 
channel each. So the choices of the basis is not unique anymore. Any two charge (or spin) fields orthogonal to $2 
and (^^4 (or 6*^4+, and 61^4-) can be used. For example, the simplest choice would be the same as Rp and used in 
the previous section. 

We then follow the hierarchy of these subleading terms in repulsive Hubbard model. Pick up the largest one among 
them and rewrite it in terms of the new fields according to Rp and Ra- After replace (pl^ , O14 , and 9"^^ by constants, 
The largest subleading term (C42 + C42) in Eq. (|4.14p becomes: 

cos V^(j)i2 cos \/4^9i2 

= cos V¥(\/2$2p + 0?4 - <t>^23) COS V^(6I^4+ + 9' 

^ cos(V2^$^ - V^cj^g) cos V^i9^+ + 9^). 

Similarly, next largest term (0^3 + 0^3) becomes: 

cos VAit^'I'^ cos V4^9'[^ 
cosiV2^'S>P + V^r23)cosV^i9^+ -9^). (4.21) 

With the perturbations like Eq. (|4.20|) and (|4.2ip . four more boson fields <f>2' '^'23^ ' ^23^ ^nd ^33^ wfll get pinned. 
Thus, as long as the RG flow is dominated by Eq. (|4.19p . the flnal phase should be CISO and the pinned bosons are 

^2 1 '^i4~' '^23' 1 ^14^1 ^14"' ^23' ^^'^ ^23~- This pinning pattern will make the pair operator Eq. (|2.40|) decay with a power 
law and Akp density operators Eq. (|2.42p and (|2.43p exponentially decay. So it's a pairing phase with no stripes. 
Since in this phase the pinned charge or spin bosons are all (/) or 9 fields, respectively, we can use other choices for Rp 
and i?o- and the result will be the same. 

Strictly speaking, the analysis in this section is only valid in the weak coupling region. Other phases might occur 
at strong coupling or in the t — J model. In the next section, we will study a different limit that may reveal some 
strong coupling physics. 
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FIG. 3: log[|(7i(i)|] v.s. log[(i£j — ?)] plot for several typical couplings. The parameters are chosen as t — ti,± = !,[/ = 0.01 
and the hole doping is 0.135. The slope of it gives us the divergent exponent of each coupling. The solid (red) lines are the 
numerical solutions of the RG equations. The dashed lines are pure straight lines as reference with the predicted slopes 1 (pink: 
(a), Ci4 (b)), 1/2 (blue: /i°4 (c)), 15/16 (green: Ci2 (d)), and (yellow: /fa (e), cJs (f)), respectively. As one can see, the 
numerical solution of C23 (f) doesn't agree with its predicted exponent. The number of couplings, whose slopes don't agree 
with the stability analysis, depends on the initial conditions. With the initial conditions used here, there are total 11 couplings, 
predicted with zero slope near the fixed ray, but have nonzero slopes in the numerical solution. However, these new term don't 
change the pinned bosons and the final phase is the same as that when these terms are irrelevant. 



V. LIMIT OF 2 DECOUPLED 2-LEG LADDERS 



One interesting limit in which it is relatively easy to understand the stripe phase is the limit of two 2-leg ladders 
weakly coupled by electron hopping and density-density interaction. Essentially this limit was discussed in Ref. [2]| . 
sub-section (VII-B-1) in the context of a 2-dimensional array of 2-leg ladders. As we discuss below, the low energy 
effective Hamiltonian in this limit is the same one describing the 2-lcg hosonic ladder which was discussed in Ref. |27| 
based on the "bosonization" for ID bosons [4lj. See Fig. ^ for illustration. 

This limit corresponds to t2.i_ and V2.i_ very small, is the following Hamiltonian: 



a=± 



and 
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(5.1) 



(5.2) 
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We set = ^3^^ = t for simplicity, but this is not essential. Thus we may begin by considering the behavior of 2 
decoupled 2-leg Hubbard ladders. Over a wide range of parameters, the 2-leg Hubbbard ladder is expected to be in 
a CISO phase |16ll26j . Introducing band bosons, i^^ly, 02^^ for the upper 2-lcg ladder (on legs 1 and 2 and v = p ot a) 
and then changing variables to 



■'lu 



± 



2, 



(5.3) 



this phase is expected to have 4P_p and pinned. The lowest Friedel oscillation wave-vector is 4fci?, corresponding 
to the 2-lcg version of stripes, namely equally spaced pairs of holes (1 on each leg) forming a "quasi charge density 
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FIG. 4: In the limit that t2^± and V2,± are much smaller than the minimum gap of the bosons, A, each 2-leg ladder is well- 
described by the CISO phase, which has pairing and Akp density oscillation. The direct electron hopping t2,± becomes an 
irrelevant process yet pair hopping, t' , generated by the higher order process will appear and Akp — 4fcF component, V' , is 
the lowest order relevant term in the interaction V2,_l. The phase is determined by the competition between t' and V' . If t' 
dominates, the system has pairing (boson superfluid) and Skp (47rpo) density oscillation in fermion (boson) language, where po 
is the average boson density. If V' dominates, the system has bipairing (boson pair superfluid) and ikp {2npo) density. 

wave" near the boundary. The 4>^p boson is. of course, gaplcss. Let A be the minimum gap for the other three bosons. 
In the limit U << t, we expect A oc texp[— const x t/U]. For U > t we expect A to be 0{t) or larger. Of course the 
lower 2 legs have a gapless boson 0^^. 

We now turn on small i2J_ and V2j_, coupling together the two 2-leg ladders. Both these interactions involve duals 
of pinned bosons, ^'^p^ and (f>±(f' and hence are ultimately irrelevant. On the other hand, a pair-hopping term, 
with amplitude t' (x t\ ^ and an interaction between the Akp density operators in the two 2-leg ladders, of strength 
V oc V2.± are generated perturbatively. [This linear dependence of V on V2,± follows since, (by analogy with the 
calculation in Appendix B), the ^kp term in the density operators for each 2-leg ladder is 0{U). Here we disagree 
slightly with Ref. [2l| which finds this interaction to be oc (V2,±)^.] Neither of these interactions involves the dual of 
any pinned bosons. For sufficiently weak t2.A. and V2,_l we may analyze the low energy theory by simply replacing all 
pinned bosons from the 2 2-leg ladders by their expectation values and writing an effective Hamiltonian for (p^^ and 
ip+p- The effective Hamiltonian describes the physics at energy scales ^ A and for it to be valid the energy scales 
characterizing the pair hopping and Akp — 4fc_F density interactions must also be <C A. It is now convenient to change 
boson variables to: 

<^± ^ [^Ip ± <^^:p]/V2, (5.4) 

since the pair hopping and Akp-'^kp density interactions only involve ^_ and its dual, 0^. The effective Hamiltonian 
at energy scales ^ A can be written: 

+t'cos{V2^(t)^) + V'cos{Vs^e-)}. (5.5) 

Here, to lowest order in t2.± and V2.±, = is simply the Luttinger parameter of the bosons on the upper 
and lower 2-leg ladders and likewise f+ = w_ is the corresponding velocity. The Luttinger parameter of the 2-leg 
ladder is expected to approach 1 at half-filling and to decrease as the density moves away from half-filling. 

We observe that the connection between (j)± and the fields in Eq. (|2.16p is not so straightforward, even in the 
limit t2,_L 0. Ignoring for the moment all interactions, when t2.i^ is strictly zero the bands come in two identical 
pairs, one member of each pair from the upper 2 legs and one from the lower 2 legs. However, as soon as t2,_L 7^ 0, 
these bands are mixed. This band-mixing is not taken into account in the present approach in which bosonization 
fields are introduced separately for the bands on the upper 2 and lower 2 legs. The present approach should be the 
correct one in the limit considered of very small t2,i, but it is non-trivial to connect the results with those obtained 
from the standard weak coupling approach. We note that a very analogous situation occurs even in the much simpler 
and well-studied 2-leg spinless fermion model. If the inter-chain hopping is sufficiently weak one normally bosonizes 
the fermions on each leg, whereas in the weak coupling limit, bosons are introduced for each band. Because of the 
exponentials entering the bosonization formulas the relationship between the "leg boson" and "band bosons" is very 
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non-linear. Nonetheless, it appears that the same phase diagram can be obtained using either approach. This can be 
seen by comparing the various features of the phases obtained using either method (26[ . 

It is not so straightforward to estimate the conditions on t2,j_ and V2,_l for this effective Hamiltonian to be valid. 
Fortunately, this is not important for our purposes. Normalizing the operators in Eq. (j5.5p so that: 

^ gi%/2¥</._(x)g-i\/2¥0_(y) ^ 1 



<g^V^e_(.)g-^v^e_(.y) > ^ ^ (5.6) 

\x - y\'^^- 

we see that t' has a scaling dimension of (energy)^^^/^^^-^ and V' has a scaling dimension of (energy)^^^-^- (after 
setting = 1). These energies scales must be much less than the cut-off scale. A, i.e. 

V' < a2-2^-. (5.7) 

Here we assume 1/4 < < 1, which is certainly true near half-filling. As mentioned above, t' (x t"^ and V' cx V2.±- 
A more complete estimate of these parameters is more difficult to make and could be quite different depending on 
whether the 2-leg ladders are in the weak or strong coupling domain. 

The phase diagram of the model in Eq. (|5.5p has been discussed in Ref. [2lj in the context of a 2D array of 2-leg 
ladders. Precisely the same model also arises from a treatment of a 2-leg ladder of spinless bosons in Ref. [27|. The 
boson annihilation operators on the upper and lower legs are represented as: 



^-^/^ cx er'^'f'li'' +..., (5.8) 



and the boson density operators as: 



1 



*^/-^^(x)4'^/-^(x) (xnb + -^d.^e^'J' + constant x {exp[i27rn,,a; + iy/A^O^'h + h.c.} + .... (5.9) 
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Here rif, is the density of bosons on each leg. Of course, it is hardly surprising that this low energy Hamiltonian 
describes a 2-leg bosonic ladder; in our low energy approximation, the fermionic degrees of freedom on each 2-leg 
ladder have been discarded keeping only the spinless pairs, corresponding to bosons, t' represents (single) boson 
hopping between the chains and V represents inter-chain boson back-scattering. It follows from Eq. (|5.7p that both 
t' and V' are relevant for 1/4 < < 1 (and in general at least one of them is relevant for all A'_). Thus one of 
(j)- and 9- boson is always gapped, yielding a CISO phase. There are two possible phases in which either or 9- 
is pinned [23|. These two phases have evident physical interpretations in the various underlying models from which 
Hf.ff arises. In the 2-leg boson model, the phase in which t' is relevant and </>_ is pinned corresponds to a standard 
ID superfluid phase in which the boson creation operator has a power-law decaying correlation function but the term 
in the boson density operator oscillating at wave-vector 27rn6 decays exponentially. On the other hand, the phase 
in which 9- is pinned corresponds to a boson pairing phase. Now the boson creation operator has an exponentially 
decaying correlation function. There is a corresponding gap to create a single boson. The 2-boson creation operator 
\E'(a;)\E'(a;) has a power law decaying correlation function as does the term in the boson density operator oscillating at 
wave- vector 27rrib. In the 2D array of 2-leg fermionic ladders, discussed in Ref. [2l|, the phase in which is pinned 
is a conventional 2D superconducting phase and the one phase in which 0_ is pinned is an incommensurate charge 
density wave phase. (The power-law decay in the single or double 2-leg ladder system is expected to become true 
long range order in the 2D system.) The physical interpretation of these phases in our model of 2 weakly coupled 
2-leg (fermionic) ladders is now also clear. The phase in which 0_ is pinned is a conventional pairing phase. Note 
that the density of bosons (average number of bosons per site in the 2-leg bosonic ladder) should be identified with 
the density of electrons (average number of electrons per site in the 4-leg ladder) . This follows since there are half as 
many sites per unit length in the bosonic 2-leg ladder as in the fermionic 4-leg ladder. Equivalently, we may identify 
the boson density with the fermion density measured from half-filling 

nb^5 = l-n. (5.10) 

Thus we see that there are no density oscillations at inS or equivalently Akp in the pairing phase. The phase in which 
0_ is pinned is a bipairing phase with stripes. 

Which phase occurs depends on K- and also the relative size oft' and V' . When K- is in the range 1/4 < K- < 1, 
where both t' and V' are relevant, we may estimate the phase boundary by the condition that the corresponding 
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energy scales, determined from Eq. (|5.7p are equal. Thus we expect the stripe phase to occur, for this range of if_ 
where: 

V > (t')[2-2^-]/[2-i/(2/f-)]_ ^5 21) 

Decreasing K- favors the stripe phase. Noting that we expect K- to decrease from 1 as we dope away from half-filling, 
it is natural to expect that, for large enough V /t', the pairing phase will occur close to half-filling and the stripe 
phase at larger doping, which is consistent with the DMRG result 0. 

At this point it is appropriate to point out that much of the same physics was discussed in two other earlier papers 
[l^l2^ . In Ref. [3 an effective 2-leg bosonic model was also discussed as an approximation to the 4-leg fermionic 
ladder. In that case the derivation was more heuristic than what appears here. A somewhat longer range interaction 
was chosen in the 2-leg bosonic model (up to separations of 3 lattice sites along the chain direction) in order to partially 
match the numerical results on the 4-leg fermionic ladder with those on the 2-leg bosonic ladder. In particular, the 
Friedel oscillations were compared in the 2 models and shown to exhibit stripes in both models at higher doping. 
The advantage of considering 2 nearly decoupled 2-leg fermionic ladders {t2,±, V2,± small) is that we can make this 
mapping more rigorous. In Ref. [l^l the limit of very small S is studied, for fixed U. It was argued that, starting with 
half-filling, at small 5, 2 of the bands remain in a gapped state with an average filling of 1/2 while the other band 
pair is doped. At higher doping, both band pairs are doped. It was assumed that each of these doped band pairs 
yields fermion pairs. At somewhat higher doping they argue that these pairs form 4-hole clusters. Friedel oscillations 
were not discussed. Although both of these papers discuss 4-hole clusters, as does the earlier DMRG work of Ref. 
0], none of them discuss the implications of such 4-holc clusters that follow from ID field theory considerations: 
exponentially decaying pair correlations and a gap to add a single pair of holes to the system. In particular, it seems 
likely that the effective 2-leg bosonic ladder model studied in Ref. [1^ was in the boson pairing phase, with a gap to 
add a single boson and exponential decay of the boson creation operator correlation function, but this point was not 
commented on. In this regard. Figures (16) and (17) of Ref. [3l are very interesting. Fig. (16) appears to be a plot 
of E{Nii) — E{Nii — 1) versus {Ni, — 1/2) /2L where iVj, is the number of bosons in the 2-leg bosonic ladder. In a boson 
pairing phase, we expect: 

E{N,)^ ^iN, + ^{-lf^ +0{llL), (5.12) 
where 2/i is the chemical potential for boson pairs and At, is the single boson gap. Thus: 

E{Ni,) - E{Nb - 1) ^ M + (-l)'^'A, + 0(l/i). (5.13) 

Both /i and A;, will evolve smoothly with density but this zig-zag structure of E{Nii) — E{Ni, — 1) is the signal of a boson 
gap, i.e. of boson pairing. Such a zig-zag is seen for the last three points in Fig. (16), implying that -Eg -I- Eiq < 2Eg 
(for L = 23) and a corresponding boson gap at J w .2 of A;, « .05t. A zig-zag is not seen in Fig. (16) at smaller 
Nb, despite the fact that the change in Friedel oscillations to stripes appears to occur at Sc ~ .125. Possibly this is 
because the boson gap is too small relative to the finite size gap to be observable for S closer to Sc- Fig. (17) shows 
the analogous quantity for the 4-leg fermionic ladder, E{Nh) — E{Nh — 2) plotted versus [Nh ~ l)/(4i) for even Nh- 
In this case no clear zig-zag is seen, which would indicate a gap to add a single fermionic pair, up to (5 = .2. Possibly 
the problem is again that the gap is too small relative to the finite size gap. This may indicate that the heuristic 
mapping is not working in great detail since the bosonic gap appears to be significantly larger than the fermionic pair 
gap. Clearly more numerical work on both 2-leg bosonic and 4-lcg fermionic models would be interesting, either or 
larger L or for a different choice of interaction parameters, to clarify whether or not a bosonic gap (and corresponding 
fermionic pair gap) exists in the stripe phase, S > 6c- 



VI. CONCLUSION 



We have taken a number of different approaches to the 4-lcg generalized Hubbard ladder based on bosonization and 
RG. We gave general arguments about possible phases based on possible ways of pinning bosons and the finite size 
spectrum. We studied particular phases from solving the weak coupling RG equations. We determined the phases 
which occur in the limit of two weakly coupled 2-leg ladders, using the connection with a 2-leg bosonic ladder. All of 
these approaches led to the same conclusion. It is not possible to find any CI SO phases that have both stripes and 
pairing. On the other hand, it is entirely possible to find phases in which stripes coexist with bipairing. Whether 
or not 4 leg ladders, for physically reasonable and numerically accessible ranges of parameters have such a phase 
remains an open question. DMRG results have suggested a phase with stripes, but have, so far, found no evidence 
for bipairing. We can see three resolutions of this paradox. 
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• Our methods are based on certain approximations: either weak couphng, or weakly coupled pairs of 2-leg ladders. 
It is entirely possible that other phases may exist for these systems which are inaccessible to these methods. 
Possibly these phases include ones with coexisting stripes and pairing. We remark, however, that these field 
theory methods have been remarkably successful in the past at describing many types of ID strongly correlated 
systems, including, for example, 2-leg ladders 31^. It would be an important discovery that they break down 



for the 4-leg ladder. 

• Possibly these systems do not really have a stripe phase in the sense that we are using. We have given a precise 
meaning to "stripes" in the limit of a very long 4-leg ladder. We mean Friedel oscillations at a dominant wave- 
vector of 4fc^ = 27rn where n can be taken to be the hole density. Existing DMRG results certainly suggest this 
but it is possible that careful extrapolation to larger systems might not confirm this result. 

• Possibly the stripe phase apparently observed with DMRG is a bipairing phase. We remind the reader that we 
define bipairing precisely to be a phase in which correlation functions of all pair operators decay exponentially 
but correlation functions of some charge 4 operators exhibit power law decay. Furthermore, such a phase has 
a gap to add one or two particles, but no gap for four particles. The limited published DMRG results have 
suggested that the decay of the pair correlation function may be power law and have not seen a gap to add 
two particles. Possibly the correlation length for the exponential decay is too large, and the corresponding 
gap to add two particles too small, to be observed so far. Further DMRG calculations could clarify this point. 
One could either study larger systems or else change the parameters of the model in an attempt to make the 
correlation length and inverse gap smaller. In this regard, numerical work on 2-leg bosonic ladders would also 
be useful to confirm that, as expected from field theory arguments, a boson pairing phase occurs in a wide range 
of parameters. As has been emphasized before, it may be crucial to include long range Coulomb interactions to 
understand stripe phases in real materials. 

We encourage further DMRG and analytical work to decide which of these possibilities is correct. Confirming any 
of them would be an important advance. 

Assuming, for the moment, that 4 leg ladders do exhibit stripes and bipairing, we speculate on the implications for 
the 2-dimensional Hubbard model. One might think that if 2-hole clusters form on 2-leg ladders and 4-hole clusters 
form on 4-leg ladders then perhaps, extrapolating to an infinite number of legs would simply give an incommensurate 
charge density wave. Such a state is perhaps not conducive to superconductivity. Stripes have also been observed 
in 6 leg ladders In this case, it appears that the number of holes per rung is 4, rather than 6, as might have 

been expected. This is suggestive of more exotic behavior than a simple CDW, closer to ideas about fluctuating ID 
conducting wires, that have been proposed for stripe phases in 2D. Developing a field theory description of this stripe 
phase in 6-leg ladders is an important open problem. 



Acknowledgments 

We are grateful to D.J. Scalapino and S.R. White for helpful discussions. We are indebted to H.-H. Lin for informing 
us of the improved analysis of RG equations prior to his publication. We also thank R. Pereira for the interesting 
discussion on Akp density operators and S. Capponi for useful correspondences. This research is supported by NSERC 
and ClfAR. 



APPENDIX A: MORE DETAILS OF THE PROOF 



In this Appendix, wc will demonstrate how to use the conditions Eq. (|2.33|) - (|2.37|) and Eq. (|2.38p - (|2.39p to prove 
that the correlation functions of pair operators Eq. (|2.40p - (j2.4ip and Akp density operators Eq. (|2.42p - (|2.43p can not 
both decay with a power-law. 

Before we start, please note that we adopt the usual convention for set theory in mathematics. Two sets A and B 
are said to be equal, if they have the same elements. 



1. Total Charge Mode Is Gapless 

In the following discussion in this subsection, we will only consider the renormalization of kpi such that the pinned 
d'ip bosons are orthogonal to 8ip. In other words, the conditions such as 2{kpi + kpj) = 2tt, won't occur but other 

possibilities such as kpi — kpj = or kpi — kp2 + kp^ — kp4 = are allowed. Then 9'^~ and {0ip — + Osp — 04p)/2 
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may get pinned but they are still orthogonal to Qip. In this case, the gapless charge mode will be (0ip,$ip) since 
all pinned 0^^ or (p'^^ fields are orthogonal to 8ip or $1^, respectively. We can reorder the transformed basis, 0^^, so 
that (j)[p = $ip and 9[p Gip. 

Let's first consider whether the pair operator in Eq. (|2.40p and Akp density operators in Eq. (|2.42p or (|2.43p could 
both have power law decay. In order to satisfy Eq. (|2.33p . VpA 'Upn = — tt, the index a in Eq. (j2.40p must be chosen the 
same as precisely one of the indices from the set {i, j, fc, 1} in Eq. ()2.42p or (|2.43p . For example if {i, j, fc, 1} = {1, 1, 2, 3}, 
then a must be either 2 or 3. If such choice of a exists (one counter example is {i, j,k,l} = {1,1,2,2}, UpA ■ Wp„ ^ — tt 
for any a), however, it implies u^a ■ v^n 7^ between the term 6aa in Eq. (|2.40p and (<^^~ ± ^Id) (|^-42p or 

(|2.43p . Thus, the condition Eq. (|2.35p is not satisfied and the pair operator Eq. (|2.40p can't coexist with any Akp 
density operator of arbitrary A:,^}, including stripes. 

Things are less trivial for the other type of pair operators Eq. (|2.4ip . The discussion depends on the situations of 
the indices set fc, Z}. Eq. (|2.42p and (|2.43p are reduced into different operator forms depending on the different 
choices of the indices and the inner products need to be discussed separately. For example, if i = j, k = I and i ^ k, 
then Eq. (|2.42p and (j2.43p contain neither charge nor spin </> fields. One can easily find that Eq. (|2.4ip can not coexist 
with stripes. We will skip the details since they are similar to what we just showed above [36| . 



2. Gapless Field Is Not the Total Charge Mode 

Now we should address the issue about the possibility of some pinned 9[p fields which arc not orthogonal to 0ip. In 
this case, 0ip can not be a basis field and thus the gapless charge mode won't be Qip anymore. As we discussed at the 
beginning of this appendix, whether a 6[p field can be pinned or not depends on the renormalization of Fermi momenta. 
Fermi momenta are important regarding both translational symmetry and the interactions in the Hamiltonian. A 
certain 9[p field can be pinned if it's allowed by symmetry and there is a relevant interaction involving it in the 
Hamiltonian. The oscillating factors of four fermion interactions only involve the Fermi momentum combinations 
like ±(kpi ± kpj ± kpk i kpi) whereas any arbitrary combination can be considered from symmetry's point of view. 
Therefore, it will be much easier to discuss the possible pinned 6*,-^ fields, not orthogonal to 8ip, through the possible 
new interactions containing in the Hamiltonian. 

In order to check through all the possible pinning patterns efficiently, our strategy is to start with one Fermi 
momentum renormalization condition so that the new interaction involving a 9lp field not orthogonal to Qip is present. 
Next, we will assume that interaction is relevant and 9'^^, as one of the basis fields, is indeed pinned. We further 
assume the phase has pairing, that is, either Eq. (|2.40p or (|2.4ip has power-law decaying correlations. So there will be 
two types of pairing to be discussed. Each type of pairing, will require that some boson fields get pinned. If another 
9^p field needs to be pinned in order to have pairing, it's allowed but one more Fermi momentum renormalization 
condition must occur. For a phase with pairing, we will have a pinning pattern for some boson fields but the rest of 
it will still be undetermined. Then the question is if it's possible to have stripes by arbitrarily choosing the pinning 
pattern for the rest of the bosons. To answer it, we have to check the correlation function of all Akp density operators 
carefully since some of them may correspond to stripes if additional Fermi momentum renormalization conditions are 
satisfied. After this is done, we repeat this procedure but start with a new 0-^ field not orthogonal to 9ip, that is, 
another Fermi momentum renormalization condition. When all the possible Fermi momentum conditions associated 
with the interactions arc discussed, wc will have completed the proof. 

Recall that a interaction can appear in the Hamiltonian if its oscillating factor becomes a constant, that is ±{kpi ± 
kpj ± kpk ± kpi) = or 27r, here i, j, k and I are arbitrary band indices (sec the oscillating factor in Eq. (C.2) in 
Appendix C). These interactions contain the charge field like ±{9ip ± 9jp ± 9kp ± 6'/p) after bosonized. We have to 
consider the possible combination not orthogonal to the total charge mode Qip. If i,j,k and I are all different, we 
only need to consider kpi + kpj + kpk — kpi =0 or 2tt since kpi + kpj + kpk + kpi = 2?™ and kpi + kpj — kpk — kpi 
doesn't result in the field not orthogonal to the total charge mode. If two indices are the same, say i = I, we need to 
consider 2kpi + kpj + kpk = 27r and 2kpi + kpj — kpk = or 27r where i,j and k are different. If three indices are 
the same, say i ~ k ^ I, we need to consider ikpi + kpj = 27r and 2)kpi ~ kpj = or 27r where i ^ j. If two pairs of 
indices are the same, then we need to consider 2kFi + 2kpj = 2-k where i ^ j. When each of above renormalization 
condition is satisfied, there will be a new interaction containing a 9p field not orthogonal to the total charge mode. In 
each case, we also allow other Fermi momentum renormalization conditions to occur so as to change the wave vector 
of density operators or pin some other 9p field orthogonal to the total charge mode. We have to discuss the cases for 
2kpi + 2kpj = 27r and 2kpi + kpj + kpk = 2t: in details. As for other conditions, there will be a general argument. 

In the following, we only show the details for the case when 2kpi + kpj + kpk = 2t: and skip the case when 
2kpi + 2kpj = 27r. For further details please see Ref. [3^. We consider all possible four fcrmions interactions to 
determine the phases (or pinning patterns). We will show that pairing and stripes still can't coexist even if the gapless 
mode is not the total charge field. 
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a. New Interactions Due to 2kFp + kpq + kpr ~ 2n 



Now we consider different renormalization conditions on Fermi momenta, 2fci?p + kpq + kpr = 27r, which can 
also result in a pinned charge field not orthogonal to the total charge mode. Without loss of generality, we choose 
2kpi + kp2 + = 2tt. In this case, for example, the interaction like '4'Ria'^R2a'^Lia'4'L3a can be present in 
the Hamiltonian because its oscillating factor becomes a constant. If this new interaction is relevant, the charge 
field (20ip + 92p + 93p)/^/6, not orthogonal to the total charge mode, will get pinned. Assume it's really pinned 
and therefore 8ip can't be an element of the pinning basis anymore. What's the gapless mode in this case? We 
know that the field 0ip is never pinned even though the total charge field is not an element of the pinning basis. 
This implies that the combination like ^4p — 9ip is also unpinned. If it's pinned, combining with the pinned field 
(26'ip + 92p + 6'3p)/V6, will imply that 9ip is pinned. One may think that the field (04p — 9ip)/\/2 should be the 
gapless mode, yet it's not orthogonal to (26'ip + ^2p + d3p)/V^, and can't be chosen as an element of the pinning 
basis. So what gapless field can let 6ip and (6'4p — 6ip)/\/2 both be unpinned? It could be either {9ip,(j)ip) or 
{9ip,4>ip). However, any pinning basis field should be orthogonal to the pinning basis field (26iip + 92p + 6'3p)/%/6. 
Thus we conclude that {9jip,(j)ip) is the gapless mode in this case. Besides (20ip + 92p + 9^p)/^/&, the interaction 

'!/;]j-^^'(/'^2aV'Lia'0L3a also Contains another charge boson field (/)23 . Then we should assume that (l)?2^ is also pinned 
by this relevant interaction. We find that orthogonal matrix representing the basis in the charge channel should be 



Rp 



f 1 \ 

2/V6 i/Ve i/Ve 

1/V2 -1/V2 
\l/Vi -1/V3 -1/V3 J 



(Al) 



With the condition that (04p, 04p) is gapless and (26'ip + 6^2p + '^3p)/V6 and 023" are pinned, the correlation functions 
of the pair operators Eq. (j2.4ip decay exponentially. The only possible pairing operator that might have power-law 
correlations correspond to a = 4 in Eq. (|2.40p . In addition we must assume that 94a- is also pinned. Now we should 
check the Akp density operators. We have to consider four situations for the indices in Eq. (|2.42p and (|2.43p . where 
the density operators reduce to different forms. These four situations are: (1) all the indices are different; (2) two 
indices are the same but different from the other two; (3) two pairs of indices are the same; (4) three indices are the 
same but different from the last one. 

For case (1): When i,j,k and I are all different in Eq. (|2.42p and (|2.43p . Eq. (|2.40p can not coexist with both Eq. 
and since u„a ■ v^n ^ between the term 9aa in Eq. (PIIO)) and (0°^^ ± ) in Eq. ^(L^ or ^L^. 

For case (2): We assume only two indices are the same in Eq. (|2.42p and (|2.43p . Then Eq. (|2.42p and (|2.43p will 
reduce to the operators with three band indices. We just have to work out all the operator forms when two out of 
four indices are the same in Eq. (|2.42p and (|2.43p . The resultant operators have the general forms: 



g-»v^[(\/2e.p+e;+)+(,/,;-)+(\/2e.„±ej+)+(0;-)]^ ^^2) 
g-iv^[(\/2e.p+e;+)+(,/,^^-)+(ej-)+(\/2,^,„-0;+)]^ l^p^^^ 
g-Jv^[(v^e.p+e;+)+(v^0.p-</.;+)+(%/2e.<,+e;+)+(%/20.„~0;+)]^ ^^^-j 



where i, j and k are all different. In fact, Eq. (|A2p - (|A5p don't cover all the possible operators forms. However, these 
other cases correspond to simply changing the signs in front of last three parentheses in Eq. (|A2p - (|A5p . We only care 
about whether the inner products between coefficient vectors are zero or not. Therefore, we can ignore the signs. 

The correlation functions of Eq. (|A4p and (jASP will decay exponentially because (20ip + 02p + ^3p)/"\/6 is pinned. 
Due to the (/)^jr term in Eq. (|A2p and (|A3p . we can only choose {j, k} = {2, 3}. The index i could be either 1 or 4. For 

1 = 4, the correlation functions of Eq. (|A2p and (|A3p could decay with a power-law if and only if [9ip — 02p — ^^3p) / \/3 
is also pinned, which requires an extra condition on the Fermi momenta but it's possible. The point here is that the 
corresponding wave-vector is kp2 + kp^ -\- 2kp4 and it only corresponds to 27rn if kpi = kp^. However, we started 
with the condition 2kpi + kp2 + ^^^3 = 27:, which implies kpi — kp4 = 27r — 27rn ^ 0. Therefore, these are no stripes 
in that case. For i = 1, Eq. (jA2l) and (jA3p correspond to the wave-vector 2kpi + kp2 + kp^ ~ 2tt, which doesn't 
correspond to that of stripes, either. 

For case (3):Next we discuss Akp density operators where two pairs of indices are the same in Eq. (|2.42p and (|2.43p . 
Recall that we are considering the pinning pattern in which (^34'', 034") is gapless and 9^2 ■,4>m, and 9^^ (or ^40-) are 
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pinned. The Akp density operators reduce to the forms: 

g-.v^(fl^7±«^)^ (A6) 
e-^^(^&^±^r), (A7) 

(A9) 



where i ^ j. In order to have power-law decaying correlation functions, these operators should contain the gapless 
mode (6'4p, (/)4p), that is, one of {i,j} should be 4. Then Eq. (|A7p . (jA8[) and (jA9p will have exponentially decaying 
correlations since they contain the dual of the pinned boson, {26ip + 02p + ^3p)/v^ or O^a- The correlation function 
of Eq. (jA6p with = {1,4} can decay with a power-law if {9ip — 02p ~ ^3p)/'\/3 is pinned. So what's the 

interaction we need to pin the field (dip — 92p — 6'3p)/\/3? We started with the interaction 4'Bia'^jli2a''J^Lia4'L3a and 
{20ip + 02p + d3p)/V6 and 

4'23 pinned. Is it possible that (26*1^ + 02p + ^3p)/'s/6 ,^^'23 and {9ip — 92p — 93p)/\/5 
are all pinned? It's possible if there is another relevant interaction involving the the fields {9ip — 92p — 9zp)l^f^ but 
not involving 023^. This means that the 9 field in that interaction must be orthogonal to (/ij;^ but not orthogonal 
to {9\p — 92p — ^?3p)/\/3- The 9 field in that interaction also has to be orthogonal to the gapless mode {94p,(l)4p). 
Then the interaction containing ^23^, such as '4'R2a'^R2a''^L3a'>pL3a-, is what we need. This interaction can appear in 
the Hamiltonian if the condition kp2 + kp^ = tt is satisfied. However, the necessary condition that the Akp density 
operator with {i, j} = {1, 4} can correspond to stripes is {kpi, kp^} = {kp2, kp^}. This will lead to the contradiction 
that kpi + kp2 + kps + kp4 = 2tt 27rn. For the choices {i,j} 7^ {1,4} in Eq. (|A6p . the correlation functions of Eq. 
(|A6p will decay exponentially since ^23" is pinned. 

For case (4): When three indices are the same in Eq. (|2.42p and p.43[) . their correlation functions will decay 
exponentially since the condition Ua-n ■ Wo-n = is never satisfied for those operators. This is in general true and 
independent of the pinning patterns and we don't have to consider this situation for Akp density operators. 

In this subsection, we have shown that pairing and stripes can't coexist when the gapless mode is (04p, (jjip) and at 
least one Fermi momentum condition 2kpi + kp2 + kps = 2tt is satisfied. We will skip the details for the case when 
there are new interactions due to kpp + kpq = tt since they are similar to the case we just showed. 



b. New Interactions Due to Other Conditions 



There arc four more types of Fermi momentum renormalization conditions that can lead to the presence of new 
interactions that may pin the boson orthogonal to Oip. They are: 

3kp, + kpj = 27r, (AlO) 

2kp^ - kpj + kpk = or 27r, (All) 

ikp, ~ kpj = or 2tt, (A12) 

kpi + kpj + kpk — kpi = or 27r, (A13) 

where the indices «, j, fc and I are all different. Although these conditions will make some interactions non-oscillating, 
they can't be relevant since they will inevitably contain some charge boson and its dual at the same time, in other 
words, Upn ■ Vpn 7^ for these interactions. 

This can be seen easily. According to our convention to define the boson fields, we know that the relation between 
chiral fermions and charge bosons are: 

(Pffia = :::^{9,p + (f>,p + ---), (A14) 



2V2 
1 

2V2' 



IfiLxa = Tn^i^ip -(kip^ )• (A15) 



When Eq. (|A10[) is satisfied, the interaction like V'ljiaV'-LiaV'ljiQ^Lia (or i and j exchanged) can appear in the 
Hamiltonian. However, this interaction contains {—Wip + 9jp) and {—4>ip + (I'jp) which are not orthogonal to each 
other. Thus, it can not be relevant since a field and its dual can't be both pinned. 

As for Eq. (jAlip . (|A12[) and (jA13[) . the corresponding interactions will have unequal numbers of right and left chiral 
fermions. For right moving fermions tpRia, is associated with the Fermi momentum +kpi while —kpi is associated 
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with the left moving fermions tl'Lia- Eq. (jAlip . (|A12p and (|A13|) are composed of three positive and one negative 
Fermi momentum. Therefore, the corresponding four fermion interactions with zero charge must contain three right 
and one left (or three left and one right) moving fermion fields. These interactions will inevitably result in Up„ • Wp„ 7^ 
and as a consequence they can't be relevant. 

c. More Than One New Interaction 

In the above discussion, although more than one Fermi momentum renomalization condition is allowed to occur, we 
always implicitly assume that the new interaction associated with the first Fermi momentum renomalization condition 
is the most relevant one. This assumption is related to how we determine the basis field. We always choose the basis 
fields guided by the interactions. For example, if the interaction ^/']j]^„'0l{iQ^-L2ct^L2a is relevant, then will be 
pinned. Whether 6'^^ is an element of the pinning basis is the issue here. If 6'^^ is not a basis field, O^^ should be 
expressed in terms of the combination of the pinning basis fields. If the interaction V'/jiaV'|jiQ^L2aV'-L2Q; is the most 
relevant interaction, there is no reason to write in terms of other fields and we will choose as an element of 
the pinning basis fields. However, what if the interactions such as ipRia''Pma''^L2a'4'L2a and '0|jiQ,V'liiaV'i3aV'L3a both 
appear in the Hamiltonian and are equally relevant? We know that 9'^2 and 0'^^ should be pinned but how do we 
choose the pinning basis? First of all, since Oip is unpinned, then 9^^ and 9!^^ are also unpinned. Then the gapless 
mode is (04p,(/)4p) in this case. Now we can obtain the other three pinning basis fields by applying an orthogonal 
transformation to the band boson fields 9ip, 92p and 03p. Rewrite and 9^^^ in terms of the pinning basis fields and 
see which pinning basis fields should be pinned. However, there are three basis fields but only two constraints (^i^ 
and 9i^ are pinned). To get a CISO phase, we need another interaction. Is it possible to pin a 4> field in this case? If 
such charge field exists, it has to be orthogonal to 6*^^ and 9'^^ . Due to the charge conservation, any pinned charge 
(j> field also has to be orthogonal to $ip. This is impossible for the cf) field obtained from the linear combination of 
4'ipi 4>2p and 03p. So the third pinned charge field is also a 9 field. 

Now comes the question: Is it possible that pairing and stripes can coexist under the condition when (04p,(?I>4p) is 
gapless and all three pinned charge boson are 9 fields? The charge fields in the Akp density operators Eq. (jA6P - (jA7P 
arc 9 fields. The pair operator Eq. (|2.40p with the choice a ~ 4, has no other 9 dependent charge fields since the 
gapless mode (/)4p is the only charge field in Eq. (|2.40p . Then the correlation functions of Eq. (|A6[) - (|A7[) and Eq. 
(|2.40[) may decay with a power-law at the same time. Indeed, this is what happens. However, in order to pin three 9 
fields, three conditions of Fermi momenta regarding the band indices {1, 2, 3}, such as kpi + kp2 = tt, kpi + kpy, = tt 
and kp2 = kp^, or kpi + kp2 — tt, kpi = kp^ and kp2 = kp^ etc, need to be satisfied. Recall that Akp = 2Tm is 
always satisfied. It also needs one additional Fermi momentum renormalization condition so that the wave-vectors 
of Eq. (|A6p and (|A7p correspond to ikp. There are five equations to solve four unknown Fermi momenta. We find 
there is no solution for all the possible cases. Thus pairing and stripes can't coexist under these situations. 

APPENDIX B: 4fcF DENSITY OPERATORS 

The higher order components of density operator is known in ID but it's not clear what's the generalization in 
ladder systems. Here we will derive the higher Fourier modes of density operators through the process of integrating 
out the large momentum modes in the perturbative fashion. The density operator on the a^^ leg can be written: 

Using Eq. (j2.5p . we decompose na{x) into components that oscillate with various phase factors kpi ± kpj. We refer 
generically to all components that oscillate with phases ±{kpi + kpj) as "2kp" terms. Naively, these appear to be 
all components of the density operators. However, there are actually additional Akp (and higher) components. These 
arise from considering more carefully the RG transformation which leads to the low energy effective Hamiltonian. This 
transformation corresponds to integrating out, within the Feynman path integral, the "fast modes" of the fermion 
fields; i.e. all Fourier modes except for narrow bands, of width A, near each Fermi point, ±kpi. We consider in detail 
how this produces 4:kp terms in na{x), in lowest order in the Hubbard interaction, U. Consider calculating some 
Green's function involving na{x), or < na{x) > with open boundary conditions. Expanding the exponential of the 
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action, to first order in U , inside the path integral effectively adds an extra term to na{x): 

. L 4 

na{x,T) na{x,T)[l + U dr' ^ r')n(,^l(a;', t')]. 

x' = l 6=1 



(B2) 



We now expand the second term, of 0{U) in band fermions using: 

na{x,T)= ^ 5„5a,e^9-^J„(p)V,„(<Z+p), (B3) 



L 4 

x' = lb=l b,ii,i2,i3,i4,pi,P2,P3 



Here: 



^21,22,23,^4 — ^aii Sai2 ^ais Sai4 : (-^^) 



and we have suppressed the imaginary time labels t, t' which are not too important. Each fermion field, "tpia^p), 
may either be a slow mode with \p — kpi\ < A or \p + kpi\ < A or it may be a fast mode with \p ± kpi\ > A. Doing 
the functional integral over the fast modes eliminates some of the fermion fields from the correction, 6na{x), to the 
density operator, na{x), replacing them by their expectation value. To generate Akp terms in na{x) we take the case 
where four of the six fields in riaHint are slow modes and two of them are fast modes. For instance consider the case 
where: 

P = -kpi + p, i = 1 

q = kpi + kp2 + kp3 + kp4 + q 
Pi = kpi +pi, ii = 4 
P2 = -kp2 +P2, ^2 = 2 

P3 kps+ps, 13 = 3 (B6) 

where all the p and Pi obey \p\ < A and q is also small, of 0(A). Then four of the fields are slow modes but ipja{p + q) 
and '!/'ij|(pi — P2 + Ps) are fast modes. Note that we have chosen the band index to correspond to the momentum 
range for all slow modes. This would be necessary if we assume that the momentum range A around each Fermi 
momentum is smaller than the difference of Fermi momenta between different bands. The product of fast mode fields 
gets replaced by its expectation value during the RG transformation: 

< ■4^ja{p + q)i^l^-^{pi -P2 +P3) >« Sjija-^^ip + q-pi +P2 -pa)- (B7) 

From Eq. (|B6p we see that the last (5-function in Eq. (|B7|) can be written: 

S{p + q - pi + P2 - P3) ^ S{p + q - pi + p2 - Ps)- (B8) 
Thus the extra term in the density operator can be written schematically as: 

Sna{x) (X U exp[i(fcFi + kp2 + kpa + kp4)x] ^ SaiSajCj423 exp[i{-p + pi - p2 + P3)x] 

j,P.Pl,P2,P3 

ll^l^i-kpi +p)lp4'^{kp4 + Pl)lljl^{-kp2 +P2)lp3l{kF3 + Ps) 
= U exp[i{kpi + kp2 + kpj, + kp4)x\SalSa2Sa3Sai'4^'^Li^^{x)ll^R4:\{x)%l;'^^^y{x)lpB3i{x). (B9) 

Naturally, a large number of other such Akp terms are generated by choosing other momentum ranges for the slow 
and fast modes. It turns out all the terms allowed by the symmetry will be generated in the low energy continuum 
limit, which is what we expected. 
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APPENDIX C: INITIAL VALUES FOR RG EQUATIONS 



Here we explicitly give the bare coupling values in Eq. (|4.2p in terms of the interactions in Eq. (|2.ip and p.2p for 
general doped A^-lcg ladders. By using ~a ap ■ ~^ ^5 ~ "^SasSis-y — ^0/3(^75, Eq. (|4.2p can be written as 

a.p ij X 

+ + i;^)<,aV'i?.a^U^L,73 - ^/5^LoV'fl./3^U^L,a], (CI) 

where i and 7 are running from 1 to A^. Note that we have a factor of 1/2 difference from the definition of operator 
Jij in Ref. [2l| . Expand Eq. (|2.ip and (|2.2p in terms of the chiral fermions Eq. (|2.5p and we have 



a/3 ijkl X Pi 



^P,io.'4'P23C.^P3kp^PiW 

+ A^^uuS-^^^^^^^^'^^ -'^*^^^'''^^'^A^^^p^.^Akii^p^n. (C2) 

ijkl X Pi 

where = ± for fermions and with the Sj„i in Eq. (|2.3p 

TV 

^jjfci = ^ ^ Sim^jm^km^lm, 
m— 1 

^ijkl — '-'mi'^mjJm+l,k'-'m+l,l- 

The 1/2 factor for V and V± in Eq. (|2.2p will make them in the equal footing as U. Now we just have to compare 
the coefficients in Eq. (|4.2p and (|C2p for the same interaction then we can obtain the bare initial values of the RG 
interactions. Recall that /y = fji, Cij = cji and fa = 0. Following the convention in Ref. [2l| . the RG equations are 
written down for ca, Cij and fij where i < j. It will be convenient to define the following quantity for OBC: 

w-i 

Sijkl = E ^'Skl- 
ni—l 

In this basis, we write down the general form for the initial values in the RG equations: 

cf, = 2[(2-cos2fc;..)yA,,,, + y^5,,,, + C/A,,,,], (C3) 

Cf^ = 2{VAiJ,JCOs2kF,+V±Su^^ + UAu^^), (C4) 

Cy- = HVA.^ij [2 cos(fcF. - kp^ ) - cos(A:f. + kp^ )] + VI 5^^^ + UA.jij}, (C5) 



c:, = 4:[VA,jij cos(A:f. + kp^) + V±Sij^j + UAij,.j\, (C6) 



/f^. = ^{VA,j,, [2 ~ cosikp, + )] + {2S,,j, - ) + [/ } , (C7) 
/5 = A[VA,j,j co^ikp, +kF,) + V^Sujj + UA,j,,]. (C8) 



where z < j here. Eq. (|C3p - (jC8p are not the basis so that the RG potential exists. In practice, we always deal with 
the RG equations in the potential basis. Therefore, we rescale Eq. (|C3P - (|C8P into the RG potential basis (without 
tilde) by 

cf, = 4^/2(2^«,04, (C9) 
~C = ^\[\{2tiv,)cI, (CIO) 
-- ^^{2^)a1^, (CU) 



a. 
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here a is c or / and again i < j. With all above results, we have the bare initial values for the RG equations derived 
from simply taking the derivative of the RG potential in Ref. [1^. 
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